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The energy of the molecules of a liquid in their equilibrium positions and the frequency of 
vibration about those positions are expressed in terms of a parameter a, the average distance 
between the molecules. This follows from a formulation of an average potential or spherical 
symmetry from an assumed force field between pairs of molecules. The ordinary thermodynamic 
quantities are calculated from the resultant harmonic oscillator-partition function. Order- 


disorder effects are neglected. 


INTRODUCTION 


HE partition function Q for monatomic 
liquids and gases? may be defined as 


(1) 


where 
0,= f f vie J WAT Medes @ 


is called the configuration integral. 

W=potential energy of the whole system, re- 
ferred to zero for infinite separation of the mole- 
cules and dr;, drz=elements of volume for 
molecule No. 1, No. 2, etc. In the harmonic 
oscillator or cell model the molecules are assumed 
to be confined to a small volume surrounding the 
equilibrium position; consequently, a configura- 
tion integral (per molecule) g, may be defined. 
Assuming spherical symmetry in the potential 
function, 


exp (— U/kT)r'dr, (3) 
0 


1 Research Assistant (1942-1943), University Research 
Institute. : 

2J. E. Mayer and M. G. Mayer, Statistical Mechanics 
(John Wiley and Sons, New York, 1940). 


where U is the potential energy per molecule and 
r is the displacement of the molecule from its 
equilibrium position. If @ is the average dis- 
tance between molecules, the integral must con- 
verge at a/2, the radius of the cell, to have 
physical significance.* 

In a real liquid the molecules may migrate 
from one part of the system to another, which has 
the effect of increasing the entropy by an amount 
R, the gas constant. To correct the cell model for 
this effect it is customary to relate the configu- 
ration integral per molecule to the configuration 
integral of the whole system as follows: 


=Q,/N!}, (4) 


being called the communal entropy factor.* 
The potential energy U can be divided into two 
parts 
U=Uo+U,. (5) 


Uo is the potential energy of the molecules in 


their equilibrium positions and U, is the potential 


*The cell arrangement of the surrounding molecules 
would be destroyed if the central molecule passed into the 
adjoining cell. 
as i Eyring and J. Hirschfelder, J. Phys. Chem. 41, 249 
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Fic. 1, Potential between pairs from Morse function, 
Eq. (16), with constants from Table I. 


energy as a function of the displacement from 
those positions. For the harmonic oscillator 
model 


=Kr/2 or 22°vmr’, (6) 


where K is force constant and » is the frequency 
defined by the equation 


2r\m 


On substituting U in Eq. (3) 


(7) 


exp (— user) exp rar (8) 


the partition function, Eq. (1), becomes 


(—)" 


T 
Xexp (-NUy/AT)( - ) (9) 


TABLE I.* Constants for Morse equation. 


Neon 


35.5X10— erg AS 
83.2X erg AB 
3.08 
4.87 10-15 
2.06 


Argon 


u=16.2X 10-8 
v=10.3X10-" 
3.83 
16.4 10715 
1.66 


Ri in A 
A in ergs 
nin 


* Values of u and v were taken from R. A. Buckingham, Proc. Roy. 
Soc. A168, 264 (1938). 


By the usual methods this leads to the following 
thermodynamic functions: 


A=NUo—3NkT In (RT /hv) — 
H=NUo—3NkT+PV, (11) 
ln v/5V. 


To apply these formulae it is necessary to have 
U; and » as functions of the volume. 


(10) 


FORMULATION OF THE CELL POTENTIAL 


Lennard-Jones and Devonshire! have calcu- 
lated the average potential energy of a molecule 
in its cell using as the potential between pairs of 
molecules, 


g=u/R®—v/R*. (13) 


The constants u and v were determined from 
experimental data on the second virial coeffi- 
cients of several gases. 

The Morse function,® 


o= A exp [— m(R—R:) ] 


—Bexp[—n(R-—R)], (14) 


formulated to express the potential energy of a 
diatomic molecule in terms of displacement R, 
from the position of minimum potential energy 
Ri, has the same general shape as Eq. (13). 
There are five constants in this form of the 
Morse function, but the condition that R: be the 
minimum requires 


A/B=n/m. (15) 


These constants can be determined from the 
constants in the Lennard-Jones function by com- 
parison of the intercepts Ro, the minimum points 


Fic. 2. Diagram indicating variables used in the averaging 
of the potential between pairs. 


4 J. E. Lennard-Jones and A. F. Devonshire, Proc. Roy. 
Soc. A163, 53 (1937). 
5P. M. ‘Morse, Phys. Rev. 34, 57 (1929). 
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Fic. 3. The average potential of an argon molecule as a function of its displace- 
ment from the equilibrium position. The full line is from Eq. (26) and the dotted 


line from the expanded form in Eq. (27) 


Ri, and the points of inflection R2, of the two 
equations. Rather than make use of the point of 
inflection the ratio A/B will be set equal to } as 
has been done by Morse. The Morse function 
then becomes 


$=A {exp [—2n(R—R:) ] 


(16) 


It may be pointed out, however, that closer 
agreement of Eq. (14) with Eq. (13) may be 
obtained by the use of the point of inflection 
condition. 

From Eq. (13), 


(17) 
(18) 


Ro=(u/v)"*, value of R at the intercept; 


Ri = (2u/v)"/*, value of R at the minimum; 


and 
$(Ri) = —v"/4u. 


From Eq. (16) 


n=In 2/(Ri— Ro) 
and ‘ 
$(Ri)=—A. 


Substituting (17) and (18) in (20) 
n=5.65/(u/v)"®. 
Equating (19) and (21) 
A=v"/4u. 


(19) 


(20) 


(21) 


(22) 


(23) 


In Table I are given the values of A, m, and R; 
for argon and neon. 

The potential curves for the Morse function 
are shown in Fig. 1. 


Lennard-Jones and Devonshire averaged the 
potential energy of a molecule as it moves on the 
surface of a sphere, whose center is the equilib- 
rium position of the molecule in its cell. The 
average potential energy due to one neighbor is 
calculated. See Fig. 2. 


$=(1/2) f "$(R) sin 0d8. (24) 


¢ is then multiplied by the number of nearest 
neighbors C (C=12 for a face-centered cubic 
structure), and divided by two to obtain the 


potential energy per molecule, U. 
U=C¢/2. (25) 


Substitution of the Morse function into (24) 
and (25) gives 


CA [= 


2n*a 


2na+1 
( ) sinh 2nr —2n cosh 2nr| 


r 


—2 exp [—n(a—R,) ] 
na+1 
| 
r 

where a is the average distance between mole- 
cules and r is the displacement of the molecule 
from its equilibrium position. In Fig. 3 are shown 
plots (in full line) of 6/A =2U/CA from Eq. (26) 


for several values of a. 
Expanding (26) in powers of r and dropping 


) sinh ny —n cosh nr||, (26) 
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Fic. 4. Average vibration frequency as a function of 
the average distance between molecules, calculated from 
Eqs. (29) and (7). 


terms containing fourth and higher powers in 7, 


CA 
[ —2n(a—R) ] 
—2 exp [—n(a—R,) }} 
CA 
exp [—2n(a—R,) ](n’a—n) 


—exp [—n(a—R,) ](n’a—2n)}r?. (27) 
This expanded form is shown (in dotted line) in 
Fig. 3. Rice® has obtained an average potential by 
calculating the effect of a given displacement on 
each of the nearest neighbors. 

Comparison of (27) with (5) and (6) sets 


CA 
[ —2n(a—R) ] 


—2exp[—n(a—R)]} (28) 


NCA 


h2 
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2amkT \ 3/2 
o=( ) exp | ~2exp 


6rkTa 


Molar Free Volume, cc 


24 2 
Molar Volume, cc 


Fic. 5. Free volume of argon as a function of the molar 
volume from Eq. (32). 


and 


CA 
exp [ —2n(a—R,) |(n?a—n) 
—exp [—m(a—Ri) (n?a—2n)}, 


(29) 


K is now expressed in terms of a, and the average 
frequencies may be calculated by Eq. (7). In 
Fig. 4 the average frequencies are plotted against 
a. Maximum frequencies calculated from the 
Debye temperatures given by Clusius’ are 
1.710" for argon and 1.310" for neon. Kane® 
has also computed frequencies as functions of 
intermolecular distances, the treatment involving 
the heat of sublimation and the elastic constants 
of the crystal. 


PARTITION FUNCTION 


If U from Eq. (27) is substituted into (3), (4), 
and (9) the following partition function results. 


The effect on the configuration integral of 
dropping fourth and higher power terms in r was 
checked by graphical integration of the expanded 
form, Eq. (27), and the unexpanded form, Eq. 
(26). As could be predicted from an inspection of 


60. K. Rice, J. Am. Chem. Soc. 63, 3 (1941). 


{2 exp [—2n(a—R,) ](n’a—n) —exp [—n(a—R)) ](n’a—2n)} 


(30) 


Fig. 3 the deviation of the two integrals becomes 
greater, the larger the values of a. In general the 
approximation is not valid for a greater than Rj. 

7K. Clusius, Zeits. f. physik. Chemie B31, 459 (1936); 
K Clusius and L. Riccoboni, Zeits. f. physik. Chemie B38, 


81 (1937). 
8G. Kane, J. Chem. Phys. 7, 603 (1939). 
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PARTITION FUNCTION 


However, it is doubtful whether the cell model 
has any significance for large values of a; conse- 
quently, the approximation does not limit seri- 
ously the legitimate use of the model.* 


FREE VOLUME 


The partition function may be set up using 
“the smoothed potential’”’ model,® according to 
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which the molecule may move in a uniform po- 
tential field in a small volume element known as 
the free volume, 2;; 


3N/2 
) exp (— NU (31) 


On comparison with Eq. (30) 


6rkTa 


{2 exp [—2n(a—R)) ](n*a—n) —exp [—n(a—R) ](n*a —2n) 


In Fig. 5 the molar free volume, Nzv,, is plotted | 


against the molar volume of argon at 80° Abs. 
Eyring*?° and his co-workers have developed the 
concept of free volume extensively, formulating 
it in terms of measurable quantities. Eyring 
estimates the molar free volume of argon at the 
melting point from the entropy of liquid argon to 
be 0.15 cc. Using specific heats he estimates a 
value of 0.35. 


kTy, CAn 


(32) 


EQUATION OF STATE 


The equation of state can be obtained by 
substituting in (12). It is more convenient, how- 
ever, to use the formula 


5InQ 5In Q éa 
= kT: 
da 


obtaining, since a#=yV and y=v2 for a face- 
centered structure, 


P=kT- (33) 


P= (exp [—2n(a—R,) ]—exp [—n(a—R,) }) 


3a? | KT 


2\a 
PVT curves for argon are given in Fig. 6. At the 
higher temperatures these curves have minima. 
Kane’ has discussed such minima in terms of the 
melting point of the solid. Simon and Kippert" 
have obtained (6P/57T)y values for argon of 38.6 
atmospheres per degree at 76 degrees and 25.0 at 
86 degrees absolute. Equation (34) gives a value 
of 30 atmospheres per degree at low molar 
volumes and 35 atmospheres per degree at higher 
molar volumes, independent of temperature. 


HEAT OF VAPORIZATION 


The heat content of a monatomic liquid was 
formulated in Eq. (11). If this is subtracted from 


* The theory should hold only for liquids at lower tem- 
peratures and higher pressures. The partition function 
may also be used for the solid state in which case the factor 
e’ would be omitted. 

* R. H. Fowler and E. A. Guggenheim, Statistical Thermo- 
dynamics (Macmillan, New York, 1939). 

10 J. O. Hirschfelder, D. P. Stevenson, and H. Eyring, 
J. Chem. Phys. 5, 896 (1937); J. F. Kincaid and H. Eyring, 
J. Chem. Phys. 6, 620 (1938). 

" F. Simon and F. Kippert, Zeits. f. physik. Chemie 135, 
113 (1928). 


(--= [ —n(a— R,) ](6n—4n?a) — (3n (34) 


exp [—n(a—R,) ](2na—2) —(na—2) 


Pressure,atm 


24 25 
Molar Volume,cc 


Fic. 6. PVT curves for argon, from Eq. (34). 
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the heat content of the gas one obtains the heat 
of vaporization. 


AH 
(35) 


The PV term in Hi, may be neglected at ordi- 
nary pressures. For argon at 80 degrees absolute 
the heat of vaporization from Eq. (35) is 1350 
calories per mole and for neon at twenty-six 
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degrees 400 calories per mole, at a=Ro, to be 
compared with the experimental values” of 1500 
and 405 calories per mole, respectively. 

We wish to express our appreciation to Pro- 
fessors W. A. Felsing and H. Eyring for their 
interest in the progress of the work. 


122,W. M. Latimer and J. H. Hildebrand, Reference Book 
of alaaaeed Chemistry (Macmillan Company, New York, 
1940). 
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Measurements of the velocity and absorption of supersonic waves in carbon tetrafluoride 


were made over a range of pressures for each of two different frequencies. From the dispersion 
and absorption curves so obtained, the specific heat for equilibrium conditions, that portion of 
the specific heat due to intra-molecular vibrations, and the mean lifetime of a quantum of 


vibrational energy were calculated. 


INTRODUCTION 


HE velocity and absorption of supersonic 
waves in carbon tetrafluoride have been 
measured over an effective frequency range from 
105.2 ke to 2,195 kc, which represents the major 
portion of a dispersion region. The measure- 
ments were actually made at two frequencies, 
290.7 kc and 615.6 kc, while the pressure of the 
gas was varied from 3 to 3 atmospheres. The 
effective frequency at atmospheric pressure is 
then given by the ratio of frequency to pressure. 
The temperature of the various measurements 
ranged from 20.5°C to 23.5°C, but during any one 
measurement it never varied more than 0.1° and 
was known to the same degree of precision. 

A very pure sample of carbon tetrafluoride was 
furnished by Dr. J. H. Simons, the director of the 
fluorine research laboratories of The Pennsyl- 
vania State College. The supersonic chamber in 
which the measurements were made had a leak 
rate, when thoroughly degassed and evacuated, 
of 1 micron per hour or less. The glass storage 
bulb in which the gas was confined between runs 
had a much lower leak rate—a few microns per 


day. The original gas was stated to contain 
certainly less than one percent and probably less 
than one-tenth percent of impurities. The above 
stated leak rates over the total period of the 
experiments would introduce less than one one- 
hundredth of one percent of impurities. These 
statements on purity were verified by a calcula- 
tion of the molecular weight of the gas remaining 
at the close of the period of measurements. The 
calculations, based on density measurements, 
gave a value of 88.2+0.2; which is to be com- 
pared to the theoretical value of 88.01. 


EXPERIMENTAL 


The measurements were made with the Pierce 
interferometer described by Telfair and Piele- 
meier.t The absorption per wave-length was 
calculated from the logarithmic decrement of the 
plate current variation. This method has received 
ample experimental** and theoretical‘ verifi- 


1D. Telfair and W. H. Pielemeier, Rev. Sci. Inst. 13, 
122 (1942). 

2W. H. Pielemeier, Phys. Rev. 34, 1184 (1929). 

3 A. van Itterbeek and P. Mariens, Physica 4, 207 (1937). 

4H. C. Hardy, Doctorate Dissertation, The Pennsyl- 
vania State College (1941). 


SPECIFIC HEATS OF CF, 


cation. Absorption and dispersion curves were 
obtained by plotting the measured values against 
the common logarithm of the ratio of frequency 
in cycles per second to pressure in atmospheres. 
The results, corrected to 22°C, are given in 
Table I and are shown in Fig. 1. 

It will be seen that the absorption values for 
runs Nos. 4, 5, and 6 are not in agreement with 
the other values which yield a smooth sym- 
metrical curve. No explanation of this discrep- 
ancy can be given unless it be that an inhomo- 
geneity of the sound field existed which, together 
with that part of the absorption due to the 
inhomogeneity of the field,®> was a function of 
pressure. 

To obtain values for use in calculations it was 
necessary to extrapolate the velocity and the 
velocity squared curves at each end until they 
became flat. These extrapolations, shown by the 
broken lines in the figure, give the following 
values: 

Vo =178.2+0.4 meters/sec., 

Vo? =31,800+100 (meters/sec.)?, 
V..=192.2+0.3 meters/sec., 
V.,? = 36,950+100 (meters/sec.)?, 


where Vo is the velocity corresponding to very 
low frequencies and V, is the velocity corre- 
sponding to very high frequencies. As an indi- 


TABLE I. 


V 2202 V22 
(meters/sec.)? (meters/sec.) 


32,060- 10 
32,460+30 
32,850-+50 
33,450+70 
34,0004 70 
34,960-+ 200 


22° 


0.202 
0.276 
0.353 
0.415 
0.458 
0.462 
0.434 


0.360 
0.429 
0.440 
0.389 
0.290 
0.224 


179.08+0.04 
180.2 +0.1 
181.3 +0.1 
182.9 +0.2 
184.4 +0.2 
186.7 


180.9 
183.2 
187.6 
189.7 
191.4 
191.9 


32,710+70 
33,550+200 
35,190+40 
36,000+40 
36,630+40 
36,820+40 


6.209 
6.341 


In runs Nos. 1-7, f=290.7 kc. 
In runs Nos. 8-13, f =615.6 ke. 


5 E. J. Pumper, J. Phys. U.S.S.R. 1, 411 (1939). Pumper 
considers the case with no internal absorption, and ob- 
tains at constant frequency the expression for the absorp- 
tion to be u=po+pe, where uo is a constant depending on 
the inhomogeneity of the sound field and y, is directly 
proportional to 1/p. Classical formulation further shows 
Me to be proportional to f/p. 


7 


® DATA FOR 2907 ke 
* DATA FOR 615.6 KC 
Temperature 22.0°C 


44 46 48 50 52 54 56 58 60 62 64 
Los 


Fic. 1. 


cation of the smallness of the error in the 
extrapolation, the agreement between the ve- 
locity and the velocity squared values so ob- 
tained may be noted. 

The value of (f/~)m, for maximum absorption, 
is 480+10 kc/atmos.; and the location of the 
inflection point of the velocity squared curve, 
(f/P)w, is 470+10 kc/atmos. (The subscript ‘‘w”’ 
is used, following some of the literature, for 
“‘Wendepunkt.”) This latter value was obtained 
from the coordinates of V,?, where V,? is given 
by 

Vi? = ( Vo?+ Vic?) /2. 


The maximum absorption, obtained from the 
solid curve, is 0.45+-0.1 per wave-length. 


CALCULATIONS 


To obtain the specific heats from the velocities 
use was made of the equation 


V?=RT/Mg(1+Rd/C), 


349 
| | | | 
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> | | | | | 
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if | 
in 
ve 
1e 
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se 
a- 
g 7 5.022 
1e 1 5.172 
3 5.462 
n- 4 5.596 
5 5.765 
6 6.007 
12 5.329 +0.2 
13 5.505 +0.6 
ce 11 5.802 +0.1 
9 +0.1 
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ed 
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where 
g=1+2bp/RT —2ap/R°T?, 
d=1+ap/R°T?, 


and “a” and ‘‘b” are the van der Waals correction 
terms. Since no equation of state data is available 
for carbon tetrafluoride an attempt was made to 
estimate values for “‘a’”’ and ‘‘b.”” A comparison 
was made among the values of ‘‘a’”’ and ‘‘b,”’ the 
critical pressure and temperature, and the 
melting and boiling points of CH4, CCl4, SiH, 
SiF4, and SnF,. The results, which are correct 
only as to order of magnitude, are: a=0.0065 and 
b=0.0022. 

The use of these values and the values of Vo 
and V,, in the above equation yields for 22.0°C: 


Co=13.5+0.5 cal./mole, 
C,=6.0+0.2 cal./mole, 


where Cp is the equilibrium value of the specific 
heat at constant volume and C,, is the effective 
value at infinitely high frequencies. The internal 
specific heat, caused by intramolecular vibra- 
tions, is given by the difference between the 
above two specific heats. 


C;=Co— =7.540.7 cal./mole. 


The value stated above for C.. agrees well with 
C,.” =3R at infinitely low pressures. 

As a check on these values the absorption per 
wave-length may be calculated. 


=0.4740.06. 


WALTER H. BYERS 


The observed value, taken from the curve, is 
0.45+0.1. 

The above value of Cp at atmospheric pressure 
is considerably larger than the value of 11.9 
predicted by Eucken and Bertram® for Co at 
infinitely low pressures. The cause of this dis- 
crepancy is undoubtedly the values of the fre- 
quencies of vibrations which they used; they 
were not at all in agreement with later values 
given by Yost? and Bailey. 

From the smoothness of the curves and values 
of C, and C;, it is apparent that carbon tetra- 
fluoride has but one relaxation time. The mean 
value of the relaxation time at one atmosphere 
has been calculated from both the absorption and 
the dispersion curve by means of the equations 


2m(f/P)mk (Ca+R) Ce 

1 Co 

2a(f/P)w Cx 

The values so obtained are 

(by absorption curve) = (7.4+0.5) sec., 
B (by dispersion curve) = (7.8+0.5) X10-’ sec., 


giving a mean value of 


8=(7.6+0.5) X10-7 sec. 


6 A. Eucken and A. Bertram, Zeits. f. physik. Chemie 
B31, 361 (1936). 

7D. M. Yost, J. Chem. Phys. 4, 325 (1936). 

8C. R. Bailey, J. Chem. Phys. 5, 274 (1937). 
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The Structure of Liquid Carbon Tetrachloride 
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University of Missouri, Columbia, Missouri 


(Received June 1, 1943) 


Monochromatic x-ray diffraction patterns of liquid CCl, at 25°C and at —20°C have been 
analyzed to obtain radial electron density distribution curves. These curves show discrete peaks 
corresponding to distances between atoms in a molecule (1.7A and 2.9A at —20°C; 1.74A and 


2.92A at 25°C), as well as other peaks (3.6A, 4.1A, and 6.4A at 


— 20°C; 4.0A and 6.3A at 25°C) 


which are not discrete and which represent distances between atoms in different molecules. 


INTRODUCTION 


PPLICATION of the generalized Fourier 

method! to x-ray diffraction patterns of 
polyatomic substances has been made for a wide 
variety of glasses? and for a few liquids.* In a 
recent article,t this type of analysis has been 
applied to the diffraction pattern of carbon 
tetrachloride, giving as a result, the radial 
electron distribution curve for this substance. 
Previous to the appearance of this article on 
liquid carbon tetrachloride, the present writers 
had carried out an extensive series of experiments 
on the same liquid at two temperatures. In addi- 
tion to presenting, herewith, the effect of temper- 
ature on the distribution curve, certain aspects of 
these curves are compared with those of the 
recently reported distribution curve for the same 
liquid. 

EXPERIMENTAL 


Most of the diffraction patterns were taken 
with Mo Ke radiation and photographic film, but 
one excellent pattern was obtained with Ag Ka 
radiation and a Geiger counter as detector of 
the scattered radiation. In the photographic 
method, a large camera of 9.20-cm radius was 
used to insure sufficiently accurate determina- 
tions of scattering angles. A total of twelve 
useful patterns were obtained in this manner at 
25°C, as well as two fairly good patterns at 
—20°C. Six prominences appeared on all films, 
and on six films, a weak seventh peak was 
measurable. Eight of the 25°C films showed no 


1B. E. Warren, H. Krutter, and O. Morningstar, J. Am. 
Ceram. Soc. 19, 202 (1936). 

SEs. J. Biscoe, A. G. Pincus, C. S. Smith, and B. E. 
Wenn Am. Ceram. Soc. 24, 116 (1941). 

3 E.g., G. G. Harvey, J. Chem. Phys. 7, 878 (1939). 

4A, SHinstein, Phys. Rev. 63, 304 (19 43). 


clearly defined seventh peak, presumably be- 
cause of under-exposure. The pattern obtained 
with the Geiger counter and Ag Ka radiation 
gave a clear indication and satisfactory measure 
of the seventh peak. 

Several types of sample were tried. A free- 
stream sample was first used and later discarded 
when it had been shown that patterns so obtained 
were indistinguishable from those using a thin- 
walled glass capillary. Because of possible ques- 
tion of the absorption correction applied to 
cylindrical samples which are not uniformly and 
completely bathed in incident radiation, several 
of the above patterns were taken with a flat 
sample contained between thin Cellophane walls. 
Tests for chemical impurity introduced into the 
CCl, by the use of Cellophane gave no indication 
of contamination. The cylindrical samples were 
1.4 mm and 1.6 mm in diameter and the flat 
sample was 0.55 mm thick. As used in all cases of 
flat sample exposures, the normal to the sample 
face made an angle of 45° with the incident beam 
in order to enhance the weak seventh peak. 

Appropriate corrections and preparations for 
curve-fitting were made in the usual manner,! 
and the final curves for the two temperatures are 
shown in Figs. 1 and 2. From these curves, the 
st(s) function was determined and subjected to 
an analysis to determine the radial electron 
density distribution as shown in Fig. 3 for the 
two temperatures. 


DISCUSSION OF RESULTS 


Perhaps the most important departure be- 
tween the present distribution curves and the one 
shown by Eisenstein‘ is in the appearance in the 
present work of discrete peaks for the C—Cl and 
the CI-Cl distances within a molecule. In a 
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critical review of work on x-ray and electron 
diffraction by gaseous CClx, Brockway‘ concludes 
that the most probable value of the C—Cl distance 
in molecules of gaseous CCl, is 1.75A, from which 
the CI-Cl distance should be 2.87A. In the 
present work, the C—Cl distance is found to be 
about 1.74A for 25°C and 1.7A for —20°C; the 
CI-Cl distance is 2.92A for 25°C and 2.9A for 
— 20°C. These values are to be compared with 
the values of 1.85A and 2.95A obtained by 
Eisenstein. 

Assuming the known tetrahedral structure of 
the CCl, molecule in the gaseous state, the areas 
under the C-—Cl and the Cl—Cl peaks should be 
638e? and 3610e?, respectively, whereas in the 
present work, the C—Cl areas are 620e? for — 20°C 
and 750e? for 25°C; and the CI-Cl areas are 
3740e for 25°C and 3950e? for —20°C. Since the 
corresponding peaks as shown by Eisenstein‘ are 
not discrete, it is difficult unambiguously to 
assign areas to them, though reconstructed peaks 
represent roughly 775e? and 4129¢e?. 

In previous work of this sort on liquids which 
‘are not monatomic, discrete peaks have been 
obtained in almost every case, and it is to be 


5L. O. Brockway, Rev. Mod. Phys. 8, 231 (1936). 


Fic. 1. Intensity pattérn for liquid CCl, at 298°K, fully corrected and fitted to 
the independent coherent curve. 


expected that this should also be true for CCl. 
Thus, discrete peaks were observed for the 
liquids phosphorus,‘ oxygen,’ nitrogen,® chlorine,® 
sulphur,!° nitric oxide,’ nitrous oxide,’ methyl 
alcohol," and ethyl alcohol.* Only in the cases of 
methyl alcohol at 25°C and sulphur at high 
temperatures were the first distribution peaks 
not discrete for every one of these molecular 
liquids. But it has been pointed out that in these 
cases, the interatomic bonds are subject to 
rupture at high temperature. It is unlikely that 
the thermal motion of the atoms within the CCl, 
molecules could be great enough to give rise to 
distribution peaks which are not discrete. 

The exact form of the distribution: curve for 
CCl, is known to be strongly dependent upon the 
exact determination of the full form of the 
intensity curve and upon the fitting of the 
intensity curve to the independent coherent 


6C. D. Thomas and N. S. Gingrich, J. Chem. Phys. 
6, 659 (1938). 

7P. C. Sharrah and N. S. Gingrich, J. Chem. Phys. 10, 
504 (1942). 

8P. C. Sharrah, Ph.D. Thesis, University of Missouri, 
Columbia, Missouri (1942). 

*C. C. Gamertsfelder, J. Chem. Phys. 9, 450 (1941). 

10N. S. Gingrich, J. Chem. Phys. 8, 29 (1940). 

11 G,. G. Harvey, J. Chem. Phys. 6, 111 (1938). 
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Fic. 2. Intensity pattern for liquid CCl, at 253°K, fully corrected and fitted to 
the independent coherent curve. , 
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Fic. 3. Radial electron density distribution curves for 
liquid CCl, at 253°K and at 298°K. 


curve. At 25°C, the use in the present work of a 
composite curve from over a dozen usable pat- 
terns should tend to make the final intensity 
curve quite reliable. The distribution curve 
shown for this case was the result of the one 
analysis made on the intensity curve for this 
temperature. At the lower temperature, — 20°C, 
two separate fittings and analyses were made. 
Both gave two discrete peaks, at the positions 
labelled C-—Cl and CI-Cl, and the first analysis 
resulted in a small peak intermediate between 
these two. An average of the two analyses is 
presented in Fig. 3 which shows a vestige of this 
spurious peak. 

Beyond the first two peaks, the prominences in 
the distribution curve undoubtedly arise from 
distances occurring between atoms in different 
molecules. None of these peaks is discrete, indi- 
cating that these distances are not fixed. Perhaps 
the best interpretation of the origin of these 
peaks can be given in terms of the explanation” 


2C, A, Coulson and G. S. Rushbrooke, Phys. Rev. 56, 
1216 (1939). 
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for similar peaks found in the case of liquid 
phosphorus which also has a tetrahedral mole- 
cule. The peak near 6.4A represents the average 
distance between the centers of adjacent mole- 
cules, whereas the subsidiary peaks between 
about 3.5A and 5A represent frequently re- 
curring distances between nearest approaching 
atoms in different molecules. The presence of 
peaks implies slight preferred orientations, proba- 
bly due to valence saturation within the mole- 
cule, and a consequent slight repulsion between 
nearest chlorine atoms in different molecules. 
In order to use the areas under reconstructed 
isolated peaks for determining the number of 
nearest-neighbor molecules, the areas under the 
subsidiary peaks should be included, as well as 
that under the large 6.4A peak, since these areas 
arise from parts of neighboring molecules. Fur- 
thermore, it seems clear that other and probably 
less well-defined peaks at distances greater than 
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6.4A should exist, to correspond to distances be- 
tween atoms on opposite sides of adjacent mole- 
cules. Hence, the way in which the curve to 
represent the nearest neighboring molecules 
should be drawn in at distances greater than 
6.4A is not well enough determined at present to 
obtain any reliable molecular coordination num- 
ber. It is only possible to set a lower limit to this 
number. Following the contour of the distri- 
bution curves from about 3.2A up to about 6.8A, 
and then sloping off along the dotted line from 
there on, gives an area corresponding to a 
minimum number of nearest molecules. On the 
questionable use of the same type of area equa- 
tion as was used for atomic scatterers, this leads 
to a coordination number of about five from the 
curves of Fig. 3. It is believed, however, that this 
number should be considerably higher in accord 
with the above interpretation of the significance 
of peak areas. 
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Here are recorded data on the magnetic susceptibilities of 
potassium and of cesium dissolved in liquid ammonia at 
240°K and 220°K. Also included are a few data on calcium 
and barium. These solutions throughout the entire range of 
concentrations were regarded as representing an electron 
gas. Such systems would permit the distributions of the 
magnetic moments of the electrons to be followed con- 
tinuously from the degenerate region of the Fermi- 
Dirac statistics to those distributions where the quantum 
and classical statistics are indistinguishable. The general 
features of the magnetic behavior of a free electron gas 


ETALS dissolved in liquid ammonia have 
the property, primary for the present in- 
vestigation, of dissociating into electrons and 
positive metal ions.! 
In dilute solution, the conduction is electrolytic 
with the electrons as the negative ions carrying 
1 This classical work of Kraus is summarized by him in 
J. Frank. Inst. 212, 537 (1931). Also in C. A. Kraus, The 


Properties of Electrically Conducting Systems (Chemical 
Catalog Company, New York, 1922). 


(Received April 30, 1943) 


were recognized, overlaid, however, by interactions charac- 
teristic of the environment of the electrons in the solutions. 
In the light of the magnetic susceptibilities as well as of 
other properties, a description is proposed for the structure 
of the solutions with special reference to the conduction 
electrons. Barium was found to dissociate into two electrons 
per gram atom and it is inferred that calcium dissociated 
likewise, but much greater dilutions would be required for a 
satisfactory demonstration. A method of considerable sensi- 
tivity and accuracy has been developed for measuring the 
magnetic susceptibilities of substances at low temperatures. 


most of the current. The conductivity changes 
continuously with concentration and becomes 
metallic in the strong solutions. When saturated, 
a solution of sodium possesses about twice the 
conductivity of mercury at room temperature. 
Along with this behavior go changes in color 
which is a pale transparent blue in great dilution, 
deepens into a dense dark blue when more metal 
is dissolved, and acquires in concentrated solu- 
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tions a metallic bronze-like luster giving the 
solution every appearance of being a liquid metal. 
When the solvent isevaporated, the original metal 
is regained. In the case of lithium and also of 
calcium the metallic crystals which first deposit 
from the solution contain ammonia bound 
stoichiometrically, Li(NH3), and Ca(NHs3)¢. The 
former even at low temperatures dissociates into 
ammonia gas and the metal powder. 

These phenomena speak strongly for regarding 
the solutions as a liquid metal whose concentra- 
tion can be varied continuously over a great 
range. Like metals, the solutions may presumably 
be conceived as an electron gas which in better 
approximation is to be counted as subject to the 
microscopic fields in the solution. 

The magnetic measurements were undertaken 
in order to follow the distribution in spin or 
magnetic moment among the electrons as a 
function of concentration and temperature. At 
high electron concentration, we would expect 
according to Pauli’s well-known work? almost all 
the electrons to have paired-off their spins, the 
degenerate state of the Fermi-Dirac statistics, 
whereas in very dilute solutions, the elementary 
magnets would be independent of each other and 
conform to those distributions of the Fermi-Dirac 
statistics which are indistinguishable from the 
distributions of the classical statistics. To achieve 
a comparable variation in spin with the high 
electron densities prevailing in metals, it would 
be necessary according to theory to raise the 
temperature above the so-called degeneracy tem- 
perature, the order of 20,000°. By diluting the 
electrons of the metal, in our experiments with 
liquid ammonia, the degeneracy temperature 
would be greatly reduced, even considerably 
below the temperature where magnetic measure- 
ments are possible. 

The degeneracy temperature is given by 


h? ( 3N ) 
where h is Planck’s constant, m is the mass of the 


electron, k is Boltzmann’s constant, N is the 
number of electrons in the volume V. 


To= 


*W. Pauli, Zeits. f. Physik 41, 81 (1927). 
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In a previous work,’ it was shown that in its 
general features, the course of the variation of the 
susceptibilities: was consistent with this outline. 
When sodium was at a concentration 0.5, its 
atomic susceptibility was of the same order of 
magnitude as that of the bulk metal. The de- 
generacy temperature for this concentration was 
about 2000°K, ten times higher than the temper- 
ature of the substance 230°K and, hence, the 
electrons were degenerately distributed. At 


| 


Fic. 1. Apparatus for measuring magnetic susceptibilities 
of substances at low temperature. 


0.002M with the degeneracy temperature at 
about 50°K, the atomic susceptibility had in- 
creased about one hundred fold and had attained 
values expected of an electron gas according to 
classical statistics. 

These results were confirmed in the more ex- 
tensive investigation of Huster and Vogt.*® 
However, the data lacked the precision required 
for interpretation, especially in very dilute solu- 
tions where the actual forces to be measured were 


28. Freed and H. G. Thode, Nature 134, 774 (1934). 


4E. Huster and E. Vogt, Physik. Zeits. 38, 1004 (1937). 
5 E. Huster, Ann. d. Physik 33, 477 (1938). 
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Fic. 2. Apparatus for preparing solutions of 
metals in liquid ammonia. 


extremely small. Since the atomic susceptibility 
of sodium proved to have less than the theoretical 
value, the measurements had to be pushed to 
extreme dilutions in order to discover the limiting 
susceptibility. For this reason, we turned to solu- 
tions of potassium whose electrons would sup- 
posedly be ‘‘freer’’ than those of sodium; the 
conductivities at the same concentrations were 
known to be consistently higher. The suscepti- 
bilities were actually found greater and a satis- 
factory precision was obtained in the measure- 
ments. Also included in this communication are 
susceptibilities of solutions of cesium, of barium, 
and of calcium. 


EXPERIMENTAL 
Apparatus 


A refinement® of the Gouy method of meas- 
uring magnetic susceptibilities was made avail- 
able for substances at low temperatures without 
appreciable decrease in accuracy. 

A vertical tube divided into equal sections by a 
horizontal partition is suspended from the arm of 
a microbalance so that the partition comes into 
the center of the pole gap of a magnet. One-half of 
the tube is filled with the solution and the other 
with the pure solvent. Under these conditions the 
susceptibility of the dissolved substance per unit 
volume in dilute solution is given directly by 


2F 
*S. Freed and C. Kasper, Phys. Rev. 36, 1002 (1930). 


solvent — 


K= 
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where K is the specific volume susceptibility of 
the solution minus that of the solvent, F is the 
vertical force exerted on the balance arm, A is the 
cross-sectional area of the tube, H; is the field 
strength at the partition, while H, is that at the 
ends of the tube. In the more concentrated solu- 
tions, account must be taken of the fact that the 
solvent material in unit volume of the solution 
is not identical in amount with that in unit 
volume of pure solvent. This correction is usually 
small and altogether negligible in the more dilute 
solutions. 

To extend the usefulness of this procedure for 
substances at low temperatures, the method con- 
sisted in suspending the glass tube G (Fig. 1a) in 
an atmosphere of stagnant hydrogen inside a tube 
H which was kept at the low temperature by the 
surrounding refrigerant liquid. A stream of dry 
hydrogen entered the stationary tube H at the 
upper end J and flowed to the outside atmosphere 
at such a rate as to keep the air with its moisture 
from descending during the experiment and 
affecting the weight of the tube G. Baffles with 
small openings could be mounted horizontally on 
the ends of the tube K which was perforated to 
allow easy entry of the hydrogen from the tubes 
J. The refrigerant, liquid ammonia, was kept in 
the Dewar vessel F. Its cap, not shown, was 
provided with connections to a manometer and 
pump for reducing the pressure and consequently 
the temperature of the refrigerant. The tubes 7, 
T contained single junction thermocouples which 
abutted into H. The temperatures at the different 
levels of the hydrogen atmosphere as read by the 
thermocouples rarely varied by more than one- 
half degree even under unfavorable conditions. 
B was a capillary tube for the passage of air into 
the liquid ammonia to lessen the bumping of the 
liquid while it was being pumped. 

The Dewar vessel F had a circular cross section 
with exception of the region C (Fig. 1b). This 
fitted into the pole gap and was therefore designed 
to separate the pole pieces as little as possible and 
yet accommodate the tubes T and B. The cylin- 
drical tubing had been here enlarged on one side 
by means of a carbon form pressing out the glass 
while it was heated and softened without 
affecting the diameter of the circular portion. 

G contained a trap at the lower end as shown 
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(Fig. 1a). By this means the bubble of ammonia 
gas occupying the space left between the liquid 
and seal would not rise to the partition when the 
tube was inverted to its normal position for 
weighing. The liquid in G extended, then, in an 
unbroken column except for the horizontal glass 
partition. 


Calibration and Measurement 


The Gouy method for determining susceptibilities re- 
quires three measurements of force, a blank run on the 
container, a calibration with a substance whose absolute 
susceptibility is known, and finally the measurement of the 
material of interest. 

Both portions of G were filled with liquid ammonia for 
the blank run. Ideally under these conditions the force 
would be zero but since the system was not perfectly 
symmetrical small correction forces were to be expected. 
The dissymetry due to the temperature gradient was so 
small that it was ignored. The extreme variation in the 
forces was as a rule not more than +0.003 mg from day to 
day. Pure liquid ammonia served for the calibration. Its 
susceptibility had been determined within an accuracy of 
one to two percent by Huster. The relative susceptibilities 
of the solute as measured by the present method were 
practically independent of the value accepted for the 
susceptibility of liquid ammonia. The absolute value was, 
however, determined by it. 

For the calibration, the upper part of the tube G was 
filled with liquid ammonia, the lower portion was left 
evacuated and the net force noted on the balance at definite 
field strengths, i.e., definite currents in the electromagnet. 
The computations were made here as well as in all the 
measurements at 15.0 amp. with the corresponding field at 
about 15,000 gauss. The value at 15.0 was taken from a 
graph of forces against current. Although the individual 
points, in the measurements of the solutions, differed from 
the average by as much as 0.005 mg, it was found that in 
a twenty-four hour period, this average did not change by 
more than +0.002 mg. This variation together with that of 
the blank produced a'‘probable error of +0.005 mg in the 
net force. y 

The tube factor defined as the ratio of the susceptibility 
to force changed less than 1 percent over the temperature 
range —33°C to —53°C. 

At lower temperatures, the weighings were somewhat 
less accurate. When the bath was pumped to obtain lower 
temperatures, the ammonia bumped. This was almost 
entirely eliminated by aspirating air through the capillary 
tube B. To obtain a uniform temperature throughout the 
refrigerating liquid it was found necessary to heat the bath 


slightly throughout its length, a function performed with . 


bare platinum resistance wire. In this manner, the tempera- 
ture of the bath could be controlled to better than a degree 
with a temperature difference of less than one-half degree 
between the top and bottom. 


Preparation of Solutions 


Potassium and Cesium 


The solutions were prepared in an all glass vacuum line 
system adapted by H. G. Thode and R. P. Metcalf of this 
laboratory from the procedures developed by Kraus and by 
Gibson and Phipps.? The metal was introduced into D at 
the point B (Fig. 2a) by distillation. A length of capillary 
had been filled with the metal and a section of the appro- 
priate length was broken off for each concentration. The 
desired amount, within ten percent, was slowly distilled 
into D by the use of free flame. Ammonia which had been 
twice dried over potassium was condensed into D. The 
condensation was complete after there was enough solution 
in D above the upper siphon so that G could be filled and an 
aliquot portion of the remainder could be siphoned into 7. 
The solution was then stirred by creating a slight vacuum 
and so bringing about vigorous boiling. The solution was 
then successively siphoned into G, T, and S. G was sealed 
off at E and the solutions in 7, S, and D were kept for 
analysis. 


Calcium and Barium 


Because of the high melting points of these metals, 
another method was devised for introducing them into D. 


TABLE I. Susceptibilities per gramatom x4. 


Free electron 
gas 
x4 108 


Force 
(corrected) K soln—Kyy, 
(mg) X108 x4 X105 


Potassium 240°K 


Cesium 240°K 


0.00470 
0.00724 
0.00700 


Cesium 220°K 


Calcium 240°K 


0.00274 
0.00531 


Calcium 220°K 
0.00245 
0.00283 


Barium 240°K 
0.00235 


7G. E. Gibson and T. E. Phipps, J. Am. Chem. Soc. 48, 
312 (1926). 
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POTASSIUM 
240°K 
220°K 
CESIUM 
240°K 
220°K 


MOLES/LITER 


The apparatus represented in Fig. 2b was sealed onto the 
vacuum line; F was joined to F on D (Fig. 2a). After the 
system had been thoroughly pumped, flamed, and flushed 
with gaseous and liquid ammonia, a piece of the metal was 
put into B and the system highly evacuated. Ammonia 
was then condensed on the metal so that the saturated 
solution which formed ran down the spiral capillary tubing 
C and through the filter plate D, a 01 Selas clay filter. 
This solution was siphoned into D through the stopcock F. 


Analyses of Concentrations 


The volumes of the solutions contained in traps S, D, 
and T were determined by distilling the ammonia into three 
measuring cells of the type described by Kraus.® After the 
ammonia had been distilled from S, D, and T, the residue 


was dissolved in water, boiled to remove ammonia and | 


titrated with 0.01N HCl. The concentrations of the solu- 
tions in D were usually about 5 percent higher than those in 
T and S which agreed with each other within better than 1 
percent. The higher concentration in D was doubtless due 
almost entirely to some undissolved metal which had con- 
densed during the distillation considerably above the solu- 


tion line. Because of good agreement in the analyses of the" 


solutions in T and S and because siphonings were rapid, the 
concentration was taken as the average in T and S. Ina 
few instances, D was averaged in if some mishap occurred 
to one of the other. 

Previous work" has shown that the densities of the 
dilute solutions were practically identical with that of pure 
ammonia. Only in the most concentrated solution of 
potassium is the difference appreciable and the correction 
was applied. The actual density was obtained by interpo- 
lation from the data of Johnson and Meyer. 

The experimental errors were incurred almost solely 

_ through the errors in weighing; the errors in analyses were 
relatively small. 


8C. A. Kraus, J. Am. Chem. Soc. 43, 749 (1921). 

°C. A. Kraus, E. S. Carney, and W. C. Johnson, J. 
Am. Chem. Soc. 49, 2206 (1927). 

10W. C. Johnson and A. W. Meyer, J. Am. Chem. Soc. 
54, 3621 (1932). 

ul E. Huster, Ann. d. Physik 33, 477 (1938). 


Fic. 3. Atomic magnetic susceptibility as a function of concentration. 


The errors in the analyses of potassium and of cesium 
amounted to less than 1 percent, in calcium 2 percent, in 
barium 3.5 percent. 


EXPERIMENTAL RESULTS AND DISCUSSION OF 
DATA 


The data for the four metals measured are 
gathered in Table I and are also presented in the 
graph (Fig. 3). 

The following features of the data support the 
idea held at the beginning that the solutions may 
be taken to represent an electron gas. 

(1). In great dilution, the atomic susceptibility 
approximates the value which would arise from 
independent elementary magnets possessing one- 
half unit of spin. 

(2). The concentrated solutions which obvi- 
ously possess strongly metallic properties furnish 
atomic susceptibilities in rough agreement with 
the presence of an electron gas free to move 
throughout the volume. For example, Huster 
found that a saturated solution of sodium had an 
atomic susceptibility of about 80 whereas 
a free electron gas in the same volume would have 
about 40 10-°. 

(3). With increase in atomic volume, the 
atomic susceptibility increased. At 240°, the 
concentration where the steep increase sets in 
occurred at about 0.01M. The ideal theory has 
0.02M as the concentration where the degenerate 
gas begins to pass into the classical distributions 
at 240°K. With increasing temperature, Fig. 3 
indicates, the difference between the correspond- 
ing concentrations diminishes. 

It is evident that the model of a free electron 
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gas must be extended to include interaction with 
the environment. Free electrons would because of 
their motion in a magnetic field show a dia- 
magnetism” cancelling one third of the spin 
magnetism. The data establish that such a 
cancellation does not occur; the two-third value 
of the spin paramagnetism indicated by the 
lower horizontal line on the left of the graph is 
exceeded and the susceptibilities continue toward 
the full value due to spin, the upper horizontal 
line. Peierls'® has raised the question whether 
Landau’s conclusions are valid for a condensed 
system where the average frequency of collisions 
between the electrons and their environment is 
high. 

The actual susceptibilities are in general lower 
than computed" for a free electron gas (last 
column of Table I) and even more drastic de- 
parture from such a model comes to light in the 
temperature dependence of the susceptibilities 
which decrease rather than increase with de- 
creasing temperature. 

To account for such a behavior, the idea 
naturally suggests itself that electrons are re- 
moved by the positive ions to form diamagnetic 
potassium molecules 


2K++2E-2K2. 


It should be noted that the formation of potas- 
sium atoms or of molecular ions would instead 
increase the susceptibilities. Huster had assumed 
the presence of diatomic molecules in solutions of 
sodium. A computation shows, however, that the 
concentration of the molecules needed for agree- 
ment with the susceptibilities appears in conflict 
with the molecular weight of dissolved sodium as 
revealed by Kraus’!® measurements of vapor 
pressures. There are other factors too which indi- 
cate that the hypothesis of dissolved diatomic 
molecules is rather forced. The very dilute solu- 
tions and the very concentrated solutions both are 
in agreement with a system consisting only of 


electrons and separate ions. The existence of the » 


metals such as Li(NHs3)4, Ca(NHs)6, Sr(NHs)e, 
Ba(NHs)¢ containing the same number of am- 
monia molecules as completely surround the 


”L. Landau, Zeits. f. Physik 64, 629 (1930). 
3 R. Peierls, Zeits. f. Physik 80, 763 (1933). 
4 N.F. Mott, Proc. Camb. Phil. Soc. 32, 108 (1936). 
%C.A. Kraus, J. Am. Chem. Soc. 30, 1197 (1908). 
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positive ions in their salts indicates strongly that: 
the ions are separated from each other in these 
metals by ammonia. We shall take the single 
point of view that throughout the entire range of 
concentrations all these solutions consist of posi- 
tive ions and electrons in ammonia and that no 
appreciable concentration of diatomic linkage 
exists between metal atoms or ions. We shall then 
inquire into the disposition of the conduction 
electrons to which the lower susceptibilities may 
be ascribed. 

The specific volume susceptibility of a de- 
generate gas is given by!®!7 


where K is the susceptibility per unit volume, y is 
the Bohr magneton, and (dZ/dE)z=£9 is the 
number of energy levels per unit energy range at 
the top of the Fermi distribution. The energy of 
interaction of the magnetic moments of the 
electrons with the external field serves to 
uncouple the pairs of electrons in the filled cells 
and send some into empty ones so that the spins 
are in line. The effectiveness of this process 
depends upon the number of cells (or levels) 
which are available to the energy of interaction. 
Therefore, the more widely separate are the 
levels, the lower is the susceptibility. 

It has been shown by Bardeen'® that when 
resonance energies are included the resulting 
decrease in the density of levels lowers the magni- 
tude otherwise computed for the heat capacity of 
metals. Magnetic susceptibilities are affected in 
the same general way. 

This resonance binding may, we propose, be 


represented as follows: 


solvent 


Ken 


solvated 
ion 


where the ammonia near the K+ corresponds to 


WW. > Zeits. f. Physik 41, 81 (1927). 

wW. Mott and H. Jones, Theory of Properties of 
— An" Alloys (The Clarendon Press, Oxford, 1936), p 
184 

18 J. Bardeen, Phys. Rev. 50, 1098 (1936). See, e.g., F. 
Theory of Solids ‘(New York, McGraw-Hill, 
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one of the molecules of solvation and the am- 
monia structure facing it represents the solvent, 
in general, but may also on occasion stand for 
ammonia of solvation of another potassium ion. 
Structures with the electron near one of the 
hydrogen atoms may be imagined stabilized by 
the many equivalent or nearly equivalent struc- 
tures which the situation may assume. The 
mobility of the electron in electrolysis may be 
viewed as the passage of the electron from one 
hydrogen atom to another with some differenti- 
ation between the ammonia molecules of solva- 
tion and of solvent. 

Also contributing to the decrease in suscepti- 
bility and especially to its decrease with the 
temperature is the interaction of pairs of electrons 
which have some resemblance to F’ centers 
known in systems composed of a metal such as 
potassium dissolved in a crystal as potassium 
chloride. In a site where a negative ion is missing, 
two electrons can lodge with some stability.!® 
This pair is probably diamagnetic. In addition to 
the interaction having analogy with F’ centers, 
the electron pairs are further stabilized through 
resonance in the same way as the unpaired elec- 
trons. 

In our present schematization another electron 
would be attached, let us say, to one of the 
hydrogen atoms on the ammonia molecule to the 
right. The interaction of these electrons to form a 
diamagnetic pair would then account for the 
decrease in paramagnetism. Here again the 
electrical conductivity would be ascribed to the 
quantum mechanical passage through barriers 
somewhat as Farkas*® assumed in an _ over- 
simplified model. In this way there appears no 
contradiction with the molecular weights as 
derived from vapor pressure measurements. 


19N. F. Mott and R. W. Gurney, Electronic Processes in 
Ionic Crystals (The Clarendon Press, Oxford, 1940), p. 128. 
20L. Farkas, Zeits. f. physik. Chemie A161, 355 (1932). 
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With further increase in electron concentra- 
tion, the increased interaction may be viewed as 
an increase in the Fermi zero-point energy which 
operates to oppose the local trapping within 
barriers. Such a process is consistent with the 
gentle increase in susceptibilities at still higher 
concentrations which Huster found. 

The lower susceptibilities of sodium compared 
with potassium (or of calcium compared with 
barium) may be linked chiefly to the relative 
sizes of these ions. The smaller ion induces a 
greater polarization in the ammonia of solvation 
with more positive charge localized on the hydro- 
gen atoms. Relative to these hydrogen atoms, the 
exchange energies are increased with consequent 
decrease in the density of energy levels. Asso- 
ciated with this influence lowering the suscepti- 
bilities there is another aspect—the deepening of 
the energy barriers would favor an increased 
stability and concentration of the electron pairs.* 

The fact that the atomic susceptibility of 
barium is definitely greater than that arising 
from one electron per gram atom is sufficient 
evidence that the ion dissociates into two elec- 
trons per gram atom which interact of course, 
with their environment. It appears highly prob- 
able that calcium dissociates in the same way 
but, as would be expected, the interaction is 
greater. 

We wish to express here our gratitude for much 
assistance to Dr. Richard P. Metcalf and to Dr. 
Warren C. Johnson whose knowledge of the 
properties of liquid ammonia systems was always 
at our disposal. 


*In a structureless medium, one expects on purely 
electrostatic grounds that the pairing of electrons would 
also take place to a greater degree with smaller ions. The 
Debye-Hiickel ion atmosphere is here replaced by the 
electron atmosphere and with increased electrostatic 
potential from the smaller ion, there would be a greater 
probability of having more than one electron at a given 
distance from the center of this ion. 
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The infra-red spectrum of trideuteronitromethane has been studied in the range 3—-25y, and 
the Raman spectrum has been analyzed. From this data and that obtained by Wells and 
Wilson on nitromethane it has been possible to determine all the fundamental frequencies 
except that of the torsion about the C—N bond, and to carry out a normal coordinate treat- 
ment from which a potential function fitting the observed frequencies of nitromethane and 
trideuteronitromethane within 2 percent has been obtained. 


S pointed out by Wells and Wilson! the 

recent growth in industrial importance of 
nitromethane and other aliphatic nitro com- 
pounds has made their thermodynamic proper- 
ties a matter of some interest. It was in hopes 
of obtaining a potential function for the nitro 
group in nitromethane which might be trans- 
ferable to other members of this series that this 
investigation was undertaken. Such a set would 
permit a priori calculation of the nitro group 


frequencies for a statistical calculation of the 


thermodynamic properties of these compounds. 


EXPERIMENTAL 
Preparation of CD,;NO, 


The nitromethane used as a starting material 
was part of the same sample as that on which 
Wells and Wilson! performed their experiments. 
It had been specially purified and presented 
them by the Commercial Solvents Corporation. 

Trideuteronitromethane was prepared from 
this nitromethane by a slight modification of the 
procedure developed by Reitz.? His method is 
essentially simply an exchange reaction with 
DO. An acetate buffer is used to prevent side 


reactions and the exchange is carried out at ca.- 


80°C to increase the rate of reaction. 

In this particular case it appeared desirable 
to prepare 5-10 cc of material fairly high in 
deuterium content, so Reitz’s experimental 
method, which involved use of a large excess of 


* Present address: The Kellex Corporation, Jersey City, 
New Jersey. 

1A. J. Wells and E. B. Wilson, Jr., ii Chem. Phys. 9, 
314 (1940), hereinafter.referred to as XIII. 

20. Reitz, Zeits. f. physik. Chemie A176, 363 (1935). 


D.O, was not employed. Instead the exchange 
was carried out using approximately equal 
weights of DxO and CH;NO:. 

The experimental procedure used here was as 
follows. For the first exchange 5 cc of 95 percent 
D.O, 25 mg of c.p. sodium acetate and 0.01 cc 
of c.p. glacial acetic acid were added to 5 cc 
of CH;NO;, contained in a Pyrex tube. This tube 
was then stoppered (glass) and shaken on a 
combination shaker-oven for 20 hours at 75- 
80°C, after which the mixture was vacuum dis- 
tilled through a drying tube containing “Dri- 
erite’”’ and collected in a tube similar to the first 
reaction vessel. Since it was found that the dis- 
tillation removed the acetic acid as well as the 
water, in order to start a second exchange fresh 
D:O and the components of the buffer were 
added to the sample collected and the process 
above repeated. Between distillations the used 
D,O was recovered by heating the CaSO, in a 
current of dry air, collecting the water in a dry 
ice trap. 

Losses of nitromethane were rather high in 
these experiments, running nearly 3 cc per dis- 
tillation by the above method. Later work by 


.Dr. H. D. Noether* showed that considerably 


smaller losses were experienced when the dis- 
tillation is carried out at reduced pressure in 
a current of air as recommended by Reitz.? The 
nitromethane lost, together with at least part of 
the acetic acid was recovered on regenerating the 
“Drierite” so the D,O collected was badly con- 
taminated. Most, if not all, of the impurities 
could, however, be removed by distillation from 


3H. D. Noether, J. Chem. Phys. 10, 664 (1942). 
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FREQUENCY IN WAVE NUMBERS 
Fic. 1. Infra-red absorption spectrum of trideuteronitromethane. Bands drawn 


with broken lines in (a) are bands of the impurity (reference 3). 


Pressures of 


-CD;NOsz: (a) 30.5 mm; (0) 8 mm; (c) ca. 1 mm; (d) ca. 0.2 mm; (e) (heated cell, 
temperature of sample tube =75°C) ca. 320 mm. 


NazCO; (anhydrous) if the temperature was not 
raised much above 30°C. 

Two series of exchanges were run to prepare 
the CD;NO: used in these experiments. In the 
first series 5 cc of CH3;NO:2 were used as the 
starting material and the exchanges were run, 
the first with 5 g 95 percent, the second with 5 cc 
98.3 percent and the third with 5 g 99.5 percent 
D.O. The second series was started with 10 cc 
of CH;NO: and exchanges were carried out with 
(1) 5 g 95 percent, (2) 5 cc 98.3 percent, (3) 5 cc 
98.3 percent, (4) 5 cc 98.3 percent, (5) 10 g 
99.5 percent D,O. 


Deuterium Content and Purity of Product 


Calculations of the deuterium content of the 
product were carried out using the distribution 
coefficient determined by Reitz: 


(D»./Hn)/(Dw/Hw) =9.78 
at equilibrium. Here 


D=fraction of total places for hydrogen occupied 
by D, 

H=fraction of total places for hydrogen occupied 
by H, 


and the subscripts m and w refer to the nitro- 
methane and water, respectively, present in the 
reaction mixture. After making some simplifying 
assumptions which should make the calculated 


values slightly low, assuming Reitz’s constant 
holds for all deuterium contents, a value of D,, 
of 0.93 was calculated for the product of the first 
run and 0.95 for the product of the second. 

Dr. W. D. Kennedy was kind enough to do a 
mass spectrographic analysis of a sample from a 
mixture of the end products of the two runs, 
2 cc of the first and 7 cc of the second. His results 
indicated the presence of about 15 percent 
CHD.,NO:2. This corresponds to a value of 
D,=0.95, indicating the calculations and the 
assumptions on which they are based were as 
nearly correct as expected. Such a high per- 
centage of CHD.NO:z present in the sample 
might be expected to show up in the spectra. 

The sample used in the investigation of 
spectra gave no test for water with anhydrous 
CuSO,. The presence of some other impurity 
is, however, indicated by the appearance in the 
infra-red spectrum of bands which are not due 
to any of the isotopic species of CD;NO:. 
Dr. H. D. Noether® has succeeded in isolating 
this same impurity from his materials obtained 
from CD;NOz2 prepared by. the same process 
used here. For further information. on this sub- 
ject, the reader is referred to his paper. 


Infra-Red Spectrum 


As in the investigation of CHs;NOz in XIII it 
was necessary to use a heated cell to investigate 
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SPECTRA OF TRIDEUTERONITROMETHANE 


the weaker bands. A somewhat different type 
was, however, developed for use in this work 
since loss of sample was definitely undesirable. 
It was found that nitromethane had little or no 
effect on rubber when heated with it for several 
hours so cells were constructed having rubber 
gaskets between the NaCl windows and the 
25-cm Pyrex tube to take care of the different 
coefficients of expansion which had caused 
windows to crack on cooling in previous cells 
where the window had been cemented directly 
to the tube. Flanged heated copper collars 2 
inches in length were used to serve the dual 
purpose of supporting the windows and pre- 
venting the condensation of nitromethane on 
them. In order to avoid the use of excessive heat 
in maintaining the cell windows at the desired 
temperature dummy NaCl windows were sup- 
ported on the outer ends of the collars to prevent 
free circulation of air. The inner metal flanges 
to which the collars were bolted were held in 
place by the flanged ends of the glass tube pro- 
tected by rubber washers. When the entire 
assembly was bolted together the outside of the 
tube-gasket-window joint was sprayed with a 
1 : 4 solution of red Glyptal in benzene to make 
a vacuum tight seal. After drying, and baking 
for twelve hours while pumping on the cell with 
an oil pump, it appeared to be satisfactorily 
tight. No trouble was experienced, in all the times 
that the cell was used, with windows cracking. 

In operation of the heated cell, the sample was 
vacuum distilled into a trap on a side arm of the 
cell, which was then sealed off under vacuum. 
The pressure of nitromethane was regulated by 
adjusting the temperature* of an oil bath in 
which the trap was immersed, In this way the 
pressure was kept at a known value within an 
estimated probable error of about 20 mm. The 
temperature of the body of the cell, which 
naturally had to be kept hotter than the trap, 
was estimated from a calibration of heater 
current versus cell temperature obtained on the 
empty cell. Since no quantitative measurements 
were desired the 10—-15°C probable error was not 
considered important. 


* The pressure-temperature data on CH;NO2 


[see E. B. Hodge, Ind. Eng. Chem. 32, 748 (1940) ] was 
used in the absence of such information on CD;NO:2 
The difference is slight at temperatures both higher and 
lower than those used in this work. 
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Cells used in the room temperature work were 
the usual 30-cm Pyrex tubes with NaCl or 
KBr windows held in place with clear Glyptal or 
Duco which have been used in previous investiga- 
tions in this series. 

The automatic prism spectrometer used in 
these investigations has been described in I‘ of 


TABLE I. 


Raman 
Intensity 
Shift* and 
(cm~!) polarization 


3018 


2300 (e, k) 
2189 (e, k) 
2116 (e, k) 
~2065 (e, k) 


3 (broad) 
4 (P) 
3(P) 

1 (P) 


i] 
i=) 


1544 (e, k) 


1390 (e, f, 1, g) 


~1250¢ (e, k?) 
~1079 (e, k?) 
1041 (e, k) 


1090 


~1000 
944 
879 

632 


578 
563 
547 


436 
412 


895 (ze, f,g,7?, Rk) 
627 (se, i, k) 


544 (te, i?, k) 


433 (se, k) 


@ See text. 

* The letters in brackets following the frequencies of the Raman 
shifts indicate the exciting lines in the Hg spectrum from which these 
shifts were observed. Kohlrausch’s nomenclature is employed. 

** (P) indicates that the line so marked is noticeably polarized. 
Other lines are apparently depolarized. 


4H. Gershinowitz and E. B. Wilson, Jr., J. Chem. Phys. 
6, 197 (1938). 
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Fic. 2. Microphotometer traces of CD;NO. Raman 
spectra. (A) Taken with Wratten No. 2A around the 
Raman tube and 0.09 g/l of Rhodamine 5 GDN extra in 
the coolant. (B) Taken with Polaroid polarizing light 
incident on the Raman tube perpendicular to its axis. 


this series. The calibration was accomplished by 
observations on the appropriate bands in CO:, 
C:H2, C2H, and CH;Br. The CaF, prism was 
used in the range 3—9u, the NaCl prism from 8 
to 16u, and the KBr prism from 14 to 25u. 
Figure 1 shows the spectra obtained at various 
pressures, and the positions of the absorption 
maxima, together with their relative intensities 
as estimated visually from the curves obtained, 
are listed in Table I. The fact that the over- 
lapping curves taken with the NaCl and CaF, 
prisms in the 1100 cm~ region do not coincide 
is due to the much greater resolving power of the 
latter there. ; 

The low intensity of the absorption due to the 
C—D stretching vibration is immediately evi- 
dent. Considering the presence of the stronger 
bands at 2765 and 2940 cm“, one might be in- 
clined to question the reported deuterium con- 
tent if it were based only on calculations. The 
other evidence presented supports the cal- 
culations, however, and makes it appear likely 
that the higher frequencies mentioned are 
combinations. 


Raman Spectrum 


The investigation of the Raman spectrum was 
carried out using the spectrograph described by 


P. WILSON 


Gershinowitz and Wilson® and qualitative polar- 
ization data were obtained by the ‘‘Polaroid”’ 
method.® On some plates a Wratten No. 2A 
gelatin filter wrapped around the Raman tube 
was used to eliminate excitation by the 24335 (h), 
24516 (2), and 24705 (k) cm™ Hg lines and 
Rhodamine 5GDN extra was used to eliminate 
the green Hg lines as described in XIII. Cu 
comparison spectra were taken on each plate. 
Table I gives the observed shifts together with 
the exciting lines from which the shifts were ob- 
served, an estimate of the relative intensities 
and indication of whether or not the shift is 
appreciably polarized. Figure 2 shows micro- 
photometer traces of two Raman plates, one 
taken without Polaroid and the other with the 
incident light polarized perpendicular to the 
axis of the Raman tube. 

In this spectrum the 2200 and 2300 cm fre- 
quencies show up quite strongly as expected. 
The 3018 cm~ shift observed on the stronger 
plates is probably attributable to the CHD2,NO, 
known to be present, although it is possible that 
it is due to the impurity.* Absorption bands of 
the impurity were observed at 3000, 1760, and 
1240 cm and a weak 1250 cm~ shift was ob- 
served in this spectrum. The observed shift is 
probably not a combination of CD;NO, fre- 
quencies since resonance with a fundamental 
appears unlikely. 

Other than the possibility of a very weak 
polarized 980 cm™ shift, the remaining lines 
observed are readily assigned to CD3;NOx. 


THEORETICAL 
Assignment of Fundamentals 


As in XIII nitromethane is treated as having 
the symmetry C2,, using the theory developed by 
Crawford and Wilson.’ The selection rules of this 
group permit the appearance in both Raman and 
infra-red of all the fundamentals except that 
due to torsion about the C—N bond which is 
inactive. Since they also permit the appearance 
of most of the combinations they are of little 
help in making an assignment. 


5 H. Gershinowitz and E. B. Wilson, Jr., J. Chem. Phys. 
6, 247 (1938). 
ase.” Edsall and E. B. Wilson, Jr., J. Chem. Phys. 6, 124 
7B. L. Crawford and E. B. Wilson, Jr., J. Chem. Phys. 9, 
323 (1941). 
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TABLE II. 


Representa- 
tion Type of vibration 


Assign 
CHsNOz! 


Frequency (cm7) 
CD;NO2 


Calculated 


Product ratio?# 
CHsNO2 CD;:NO2 Oo 


2965 


C—H stretching 
1413 


CH; bending 

A, N—O stretching 
C—N stretching 
NO: bending 


Az Torsion about C—N 


C—H stretching 
N—O stretching 
B, CH; bending 
CH; rocking 
NO: rocking 


C—H stretching 
B CH; bending 
CH; rocking 
NO: rocking 


2965 
1426 
1363 

920 


1A, J. Wells and E. B. Wilson, Jr., J. Chem. Phys. 9, 314 (1940). 
20. Redlich, Zeits. f. physik. Chemie B28, 371 (1935) 


3 W. R. Angus, C. R. Bailey, J. B. Hale, C. K. Ingold, A. H. Leckie, C. G. Raisin, J. W. Thompson, and C. L. Wilson, J. Chem. Soc. 971 (1936). 


The assignment of frequencies to the vibra- 
tions having the symmetry A, is readily accom- 
plished by considering the polarization data on 
the stronger Raman lines. The 2189, 2116, and 
2065 cm shifts must be the symmetrical C—D 


TABLE III. 


Observed 
~1250 
~2065 
2116 
~1000 
~1276 
1310 


1527 
1545 


Assignment 
2X 627 =1254? 
2X 1041 = 2082 
2X 1079 = 2158 
560+424 =984? 
2X 632 = 1264 
879+ 424 = 1303 


944+ 632 = 1576? 
879+ 632 =1511 
~879+632=1511 


1051+ 632 = 1683 
1080+ 632 =1712 


2X 879 (2A,, 2B, and 
A,+B,)=1758 


879 +944 = 1823 


1571+424=1995 
1051+944 = 1995 


2X 1051 =2102 
2X 1090 = 2180 


1571+ 879=2450 
1384+ 1090 = 2474 


2X 1384 =2768 
1571+ 1384 = 2955 


Infra-red 


~1696 


1763 
1785 


1820 
1992 


2064 
2107 
2442 


2765 
2940 


See text. 


stretching resonating with 21079 and 21041. 
By comparison with other molecules the 1079 
cm~ shift must be the symmetrical bending of 
the CD; group, but here, in contrast to the other 
Raman shifts, the value of the gas phase infra- 
red frequency corresponding is not immediately 
evident for the 1050 cm“ is not visibly split into 
two components as it is in the Raman spectrum. 
It is possible that the parallel and perpendicular 
bending frequencies show markedly different 
shifts in going from the liquid to the gas phase, 
but since the 2000-2200 cm region remains 
practically unchanged on this transition (three 
lines at approximately 2190, 2110, and 2065 cm= 
in both spectra) it seems more probable that a 
difference still exists between the gas phase 
frequencies of the two types of vibration, and the 
1090 cm band observed using the fluorite 
prism has been somewhat arbitrarily chosen as 
the symmetrical fundamental. The assignment 
of the other symmetrical fundamentals requires 
no explanation. 

Most of the other fundamentals may be 
assigned frequencies on the basis of comparisons 
with other compounds containing analogous 
groups or from consideration of the mechanics 
of the vibrations. The broad line at 2290 cm™ 
must be the unsymmetrical C—D stretching 
frequencies, since it is the only strong unpolarized 
Raman shift in that region. 1571 cm is asso- 
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TABLE IV. 


Constant Value 


Coordinates involved* 


4.96 X 10° dynes/cm 


0.04 
0.46 
0.01 
0.665 
—0.035 
—0.35 
4.74 
9.72 
2.54 
0.636 


2.263 X — 


0.404 
1.264 


dyn 


—0.353 X 10-3 
radian 


0.749. 


Ar; in CHs 

Ar; in CH; 

ro(Aa) in CH; 

ro(Aa) -7o(AB) in CH; 

ro(AB) in H3C— NO» 

1o(AB;) -ro(ABi+1) in 
ro(Aa:) AR, 


dyne-cm 
(radians) 


Ay: Ay 
AR;(Ay;) = AR;(Ay; +i) 
AR(Ayi) 


* The coordinates involved in calculations on the methyl group are those used in previous work (see reference 9 and 16). In the table above the 
nomenclature for the remainder of the molecule is as follows: R=C —N bond, R=N —O bond, y =C —N —O angle, 5=angle between the C —N 


bond and the plane of the NOz group. 
t See reference 16. 


ciated with the NO: group and obviously is the 
unsymmetrical stretching frequency. We have 
assigned both the B; and Bz CDs; bending fre- 
quencies the value 1051 cm™ since the 1000 
cm line, if present, is almost certainly too low, 
and the two are apparently not resolved in either 
spectrum. 

The bands at 424 and 560 cm may be 
assigned to the NO: rocking motions in B; and 
Bz, respectively, on the basis of the same moment 
of inertia of the NO» group arguments used in 
XIII. This assignment is checked by the ob- 
served band envelopes. From the moments of 
inertia and the envelopes published by Badger 
and Zumwalt,® using their method, it was found 
that in both compounds the B, bands should be 
doublets while the B, bands should have ‘Q”’ 
branches. In the latter case the “Q” branch 
should be considerably weaker in CD;NO, than 
in CHsNOsz. A comparison of the spectra of the 
two compounds shows this to be the case. 

The least definitely established assignment of 
frequencies for this molecule is that for the two 


8R. M. Badger and L. R. Zumwalt, J. Chem. Phys. 6, 
711 (1938). 


CD; rocking motions. Qualitatively one might 
expect that the frequencies should be split due 
to the difference in moments of inertia of the 
NOz group about the two axes as discussed in 
assigning the NO» rocking motions, and by the 
same argument used in that case the B, rocking 
should have a lower frequency than the By». It 
also might be predicted that both should be 
lower than 970 the value found for 
by Stitt,® since the NO» group has a considerably 
greater moment in CD;NO; than the CD; group 
has in C2Dgs. If, then, the doubtful infra-red 
band at 1000 cm“ is eliminated on these grounds, 
and the 944 cm~ infra-red band is assigned to the 
B, rocking motion, the other must be lower. 
No band or Raman shift is separated, but the 
lack of a definite envelope for the 879 cm™ band 
indicates that it may actually be double. There 
seems to be no logical alternative to this selection 
on the basis of the evidence available at present 
and it is at least partially justified by the fit 
obtained with it when the product rule is applied. 
(See Table IT.) 


°F. Stitt, J. Chem. Phys. 7, 297 (1939). 
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Table III gives a tentative assignment of the 
combinations appearing in the Raman and infra- 
red spectra, on the basis of the above assignment 
of fundamentals. The doubtful band at 1000 
cm! may be a combination if C2, selection rules 
are violated, but it seems more likely that it is a 
low frequency branch of the 1051 cm band, or 
it is due to an impurity. 

Except for the slim possibility that the 1000 
cm~! band may be a combination involving it, 
there is no evidence as to the value of the 
torsional frequency in either of the spectra of 


The Vibrational Potential Function 


In the calculations of moments of inertia to 
obtain band envelopes as well as those discussed 
in this section, the dimensions determined by 
Brockway, Beach, and Pauling’® for CH3;NO2 
were used. These are 


C—N=1.46+0.02A; N—O=1.21+0.02A 
O—N-—O angle=127+3°. 


As in XIII we assume that C—H=1.08A and 
that all the methyl group angles are tetrahedral. 

We also make the usual assumption that 
neither the dimensions nor the potential function 
of the molecule is changed by substitution of D 
for H. The frequency shifts are then entirely 
due to mass differences. There is a slight differ- 
ence in the anharmonicities of the C—H and 
C—D stretching vibrations and the deviation of 
the observed from the theoretical product rule 
ratios": is in the expected direction® and is 
believed to be due to this cause. In the case of 
the B, frequencies, where the observed deviation 
is smaller than expected from comparison with 
other cases, the difficulty probably lies in the 
879 cm rocking frequency which is hidden in 
the strong C—N stretching band. To obtain a 
better fit with the expected ratio the rocking 
frequency would have to be raised slightly, which 
would produce somewhat better agreement with 
the observed combination bands as well. 


10L. O. Brockway, J. Y. Beach, and L. Pauling, J. Am. 
Chem. Soc. 55, 2693 (1935). 

uO, Redlich, Zeits. f. physik. Chemie B28, 371 (1935). 

SW. R. Angus, C. R. Bailey, J. B. Hale, ce Ingold, 
A. H. Leckie, C. G. Raisin, ji . Thompson, and C. L. 
Wilson, J. Chem. Soc. 971 (1936). 

30. Redlich, J. Chem. Phys. 7, 856 (1939). 


The normal coordinate treatment from which 
the potential function given in Table IV was 
obtained was’ carried out using methods pre- 
viously described. The secular equation was 
set up using the bond vector method and the high 
frequencies then split out.“ It was not found 
possible to factor out the N—O stretches, but 
the C—H and C—D stretches were treated 
separately, the success of the treatment being 
indicated by the fit obtained with the observed 
frequencies.* The force constants given in 
Table IV for the C—H stretch were determined 
from the frequencies observed in CH;NOz since 
work by Crawford and Brinkley'® has shown that 
in CH compounds factoring out the C—H 
stretch introduces an error of only about 0.1 
percent if no other frequencies in the molecule 
under consideration are above 2000 cm. It 
might be expected that the approximation would 
not hold so well in a case such as CD;NO2 where 
the frequency factored out is around 2200 cm~! 
and other frequencies in the molecule are as high 
as 1500 cm~. As shown in Table II, however, 
the error is only what would be expected from 
the anharmonicity of the C—H stretches. In 
comparison with the anharmonicity effect the 
error introduced by factoring out the C—D 
stretches appears to be negligible. 

When the potential function given in Table IV 
is factored as indicated it reduces to the same 
form as that given in XIII for the B,; and By 
representations. As will be noted the fit obtained 
with the values of the constants reported here is 
not quite as good as that given in the previous 
paper, especially in B,.** A better fit could un- 
doubtedly be obtained by introducing more inter- 
action terms, in particular one between the N —O 
stretch and the CH; bend, but since the agree- 
ment is probably good enough for most purposes 
and since the potential function in the two un- 
symmetrical representations treated has the 


4 E, B. Wilson, Jr., J. Chem. Phys. (a) 7, 1047 (1939); 
(6) 9, 76 (1941). 
16M. A. Er’ yashevich, (a2) Comptes rendus Acad. Sci., 
U.S.S.R. 28, 604 (1940) ; (b) J. Phys. Chem. (U.S.S.R.) 
14, 1381 (1940); (c) Ibid. 15, 831 (1941). 
* See Table II. 
6B. L. Crawford, Jr., and S. R. Brinkley, he J. Chem. 
Phys. 9, 69 (1941). 
* This is due in part to a slight numerical error made in 
rage Ny up the secular equation used in the work reported 
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advantage of almost surprising simplicity, no 
attempt has been made to improve things except 
by changing the numerical values of the force 
constants involved. 

In the A; representation the most general 
potential function, after factoring out the C—H 
stretches, would contain ten parameters. From 
the other representations it would appear that 
the interactions between motions of the NO» 
group and the CH; group should be small. 
For this reason the two parameters correspond- 
ing to such interactions were omitted. Inter- 
actions between the C—N stretch and the other 
motions are probably present.'® The other inter- 
action appearing in this representation is one 
between the N—O stretching and the NOs: 
bending. In B, the corresponding interaction is 
apparently not important. When it was omitted, 
however, preliminary investigation showed that 
the values of the other constants necessary to 
fit the observed frequencies would probably be 
so widely divergent from those expected as to 
cast doubt on the validity of the potential 
function. By including the R—y interaction it 
was possible to obtain the agreement between 
observations and calculations given in Table II 
and still keep values of the other constants 
somewhere near what was expected.* The R—+ 
interaction was chosen rather than either of the 
other two possible ones because it was felt that 
the charge on the nitro group might well lead 
to some interaction not accounted for entirely by 
“diagonal” valence type force constants. 

It will be noted that while eight parameters 
have been adjusted in this representation, there 


* By analogy with results from other compounds, it was 
expected that the C—N stretching constant, KRNO2 
=4,0-5.5 dynes/cm; also Kyy=—0.3 to —0.5 dyne/cm. 
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are actually only seven independent conditions, 
for while there are eight frequencies the assump- 
tion of identical potential functions for the two 
molecules imposes one condition on them. The 
potential function given, then, is certainly not 
unique. It does fit the facts considered in ob- 
taining it, but further tests, perhaps on other 
nitro compounds, should be made before a great 
deal of reliance is placed on it. 

It is interesting to note that if the A; repre- 
sentation is excluded eighteen frequencies have 
been fitted with a maximum deviation of 2 
percent by adjusting only five constants, taking 
the others necessary from Stitt’s work on ethane.® 
At least in these particular cases the form of the 
potential function and the numerical values of 
the constants used is quite satisfactory. 

The fact that it has proved possible, using a 
tetrahedral structure for nitromethane, to ob- 
tain a reasonable potential function which fits 
the observed frequencies of both CH;NO~2 and 
CD;NOz supports the view expressed in XIII 
that this is at least the form present in by far 
the largest amounts.'’? No frequencies appearing 
in either the Raman or the infra-red spectra of 
CD;NO: require postulation of a second form. 

In conclusion the writer would like to ack- 
nowledge his indebtedness to Professor E. Bright 
Wilson, Jr., not only for suggesting this problem 
and for his continued interest and many helpful 
discussions during the course of the investigation, 
but also for making available the original cal- 
culations on CH;NO, by himself and Dr. A. J. 
Wells. I should also like to express my apprecia- 
tion to Dr. Wells for permission to use his work 


and for his interest in this problem. 


17 See E. Pendl, A. W. Reitz, and R. Sabathy, Proc. Ind. 
Acad. Sci. 8, 508 (1938) (Raman Jubilee Volume). 
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The infra-red spectrum of cyclobutane in the range 2 to 25u has been obtained, together with 
the Raman spectrum of the same material. Analysis of these spectra, together with a normal 
coordinate treatment of the problem, has led to a tentative assignment of the fundamentals 
which explains all features of the observed spectra. Butene-1 has been identified as a product of 


the liquid phase photolysis of cyclopentanone. 


HE vibrational spectrum of cyclobutane is 

of interest on several counts. It provides an 
example of a comparatively complex molecule of 
high symmetry in which there is no internal 
rotation. Also, according to the usual model, it is 
somewhat strained and there is a possibility of 
investigating the effect of this strain on the 
valence type force constants. Recently some 
question has been raised concerning the structure 
of hydrocarbon rings, so the structure of cyclo- 
butane may require proof. There is also the 
possibility of using the spectrum obtained to 
calculate thermodynamic properties of this 
material. 

Recent work! in these laboratories on the 
photolysis of cyclic ketones made it appear 
possible that this method could be used to pre- 
pare cyclobutane in small quantities somewhat 
more simply than it can be done using the method 
devised by Willstatter? and lately modified by 
Heisig.*? This work was undertaken to test the 
practicability of the method of preparation as 
well as to investigate the spectrum of the 
product. 


EXPERIMENTAL 
Preparation of C,H; 


As indicated above, the cyclobutane used was 
prepared by the gas phase photolysis of cyclo- 
pentanone in a manner similar to that described 
by Benson and Kistiakowsky.! Certain changes 
were made in the construction of the apparatus, 


* Present address: The Kellex Corporation, Jersey City, 


ew Jersey. 

1S. W. and G. B. Kistiakowsky, J. Am. Chem. 
Soc. 64, 80 (1942). 

? R. Willstatter and J. Bruce, Ber. 40, 3979 (1907). 


5G. B. Heisig, J. Am. Chem. Soc. 63, 1698 (1941). 


however, in view of the fact that our primary 
interest was in only one of the products of 
decomposition, and as much as possible of that 
was desired. 

Since it has been shown that the relative 
amounts of decomposition products are inde- 
pendent of temperature and pressure over the 
immediately available ranges, it appeared that 
the only way to increase the rate of formation of 
cyclobutane was to increase the intensity of 
illumination. The arc‘ used in the previous work 
is the most intense source of ultraviolet radiation 
with which we are acquainted, but in those ex- 
periments the light was not focused particularly 
well. We, therefore, enclosed the arc and quartz 
reaction tube in an elliptical stainless steel re- 
flector, the arc at one focus, the tube at the other. 
An air blast had to be turned on the outside of the 
reflector to prevent its staining from the heat of 
the arc. Since the arc had to be mounted hori- 
zontally, the tube had to be similarly oriented. 

An oven was used as before to keep the leads to 
the reaction tube hot. For simplicity in con- 
struction the reaction tube was also inside the 
oven, so about half the reflector was too. Actually ~ 
no heat other than that from the arc was neces- 
sary to keep the reaction tube above the boiling 
point of cyclopentanone. 

The quartz spiral manometer was eliminated 


‘since knowledge of the true pressure in the 


reaction chamber was unimportant. The mercury 
regulator on the other side of the condenser gave 
a sufficiently accurate measure. A thermocouple 
was used in preliminary experiments to determine 
very approximately the temperature inside the 
reaction tube. 


4 Engineer’s Bulletin No. 53, General Electric Company 
(Nov. 1939). 
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The system of taking off reaction products was 
also changed. The gases, after passing through 
the reaction tube, were fed into the center of a 
vertical condenser, the cyclopentanone con- 
densing and running back into the pot while the 
products not condensable at 5-15°C collected 
above and leaked off through a mercury pressure 
regulator. Gas coming off through the regulator 
was passed through a trap in dry ice-acetone to 
remove residual cyclopentanone and any other 
high boiling substances which might be present, 
and then through a trap in liquid nitrogen to 
collect the cyclobutane and ethylene. CO was 
pumped out continuously with an oil pump. The 
photolysis and collection using this apparatus 
were completely automatic, once started. 

The cyclopentanone used as a starting material 
was a combination of samples, for one of which 
we wish to thank Dr. J. J. Leavitt, of these 
laboratories, who presented it to us, the other 
from Eastman Kodak Company. It was purified 
by distillation, first through a small Widmer 
column (b.p. 129.1-129.3°C), then through a 
3’X1” column with wire gauze (Purdue) packing 
at a reflux ratio of 30 : 1. The sample used had a 
constant boiling point of 129.2°C at 762 mm. It 
was stored in a black painted, tin foil wrapped 
flask attached to the line, and portions were 
vacuum distilled into the pot for the photolysis. 

As in the previous experiments the reaction 
system was evacuated with a Hg pump before 
starting the photolysis. Runs were made at a 
pressure of 200-300 mm and temperature of 
250°-300°C in the reaction tube. 

For larger amounts of products than those 
obtained by Benson and Kistiakowsky,' simple 

vacuum distillation is a somewhat more con- 
venient method of purification than that em- 
ployed by them. We used this, carrying out the 
distillations at the temperature of solid COs. 

Using the apparatus and conditions given 
above we were able to obtain about 0.1 cc of 
liquid cyclobutane per hour of operation at best. 
One of the factors limiting the rate of decompo- 
sition of cyclopentanone seemed to be the rate at 
which it was possible to distil it through the 
circulating system without producing bumping. 
Formation of a brown deposit, presumably a high 
polymer, on the walls of the reaction tube may also 
have cut down the effective radiation somewhat. 
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Purity of Product 


Two lots of cyclobutane were prepared, one of 
which was used for preliminary infra-red work 
and the Raman spectrum, the other for final 
infra-red work and the spectroscopic investiga- 
tion of impurities. The first had a vapor pressure 
range at —78°C of 4-8 mm. The second, after 
bromination, had a range of 4.5-7 mm at the 
same temperature (cf. 5.1 mm‘). 

Physical tests of the purity were run on the 
first batch. The tests used were the vapor pres- 
sure at various temperatures, the differential 
vapor pressure test of purity and the vapor 
density (for both identity and purity of product). 

Vapor pressure measurements at 0°C and 
— 14°C showed the vapor pressure of the sample 
to be slightly above that for pure cyclobutane.* 
This was confirmed by differential vapor pressure 
measurements by Shepherd’s method.® By this 
latter method it was estimated that 0.5 percent 
of ethylene was present, assuming that to be the 
only impurity. From measurements of the mo- 
lecular weights at various pressures® of the first 
and last halves of an approximately isothermal 
distillation, assuming the impurity to be ethylene 
and the apparent molecular weight of cyclo- 
butane to be a linear function of the pressure, a 
value of 0.35 percent ethylene was calculated. 
The measured molecular weights, when corrected 
for the presence of ethylene, plotted on a pressure 
scale and extrapolated to P=0, gave a molecular 
weight of 56.10+0.1. No evidence indicating the 
presence of either high boiling or high molecular 
weight impurities was obtained. 

Spectroscopic tests on this sample gave no 
evidence for the presence of ethylene. A band at 
728 cm— was observed but not accounted for. In 
the second sample prepared by the same method 
this line was not observed so it is almost certainly 
due to an impurity. 

In preparing the second sample fewer pre- 
cautions were taken to prevent condensation of 
cyclopentanone in the reaction chamber and as a 
result some liquid was present during the 
photolysis. The product obtained was contami- 
nated with butene-1, which was identified spec- 
troscopically by comparison with infra-red ab- 


5 M. Shepherd, Bur. Stand. J. Research 12, 185 (1934). 
6 Rayleigh, Proc. Roy. Soc. 55, 340 (1899). 
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sorption data on this compound very kindly 
supplied us by the Shell Oil Company, Incorpo- 
rated, Houston Refinery Research Laboratory.* 
While their data is only tentative the stronger 
bands seem fairly definitely established. Bands 
were observed at 1000 cm~ and 1652 cm~ in our 
impure cyclobutane (and only in that sample). 
Both correspond to strong butene-1 bands. After 
a gas phase bromination of the sample’ these 
bands disappeared, and bands at 627, 902, and 
1453 cm in cyclobutane changed shape, un- 
doubtedly due to removal of the butene-1 which 
has strong bands at 640, 910, and 1462 cm! as 
well as those mentioned above. No bands due to 
cis- or trans-butene-2,* isobutene** or methyl 
cyclopropanef were observed at any time. Since 
the bands of butene-1 were observed only in the 
sample where there was a possibility of liquid 
phase photolysis, presumably it is a product of 
that process. 

In investigating the presence of butene-1, ab- 
sorption spectra were taken of the two halves of a 
distillation of the impure sample from solid COz. 
A distinct difference in the intensities of the lines 
attributed to butene-1 was observed, as would be 
expected since butene-1 has about four times the 
vapor pressure of cyclobutane under these condi- 
tions (23.5 mm at —77.5°C*). Intensities of lines 
not corresponding to butene-1 remained constant 
within the error of measurement. 

To summarize, then, in the samples whose 
spectra are reported here there is evidence for the 
presence of ca. $ percent of ethylene. No other 
impurities were identified but an infra-red band 
observed at 728 cm™ is due to an impurity not 
present in the sample on which the greater part 
of the infra-red work was done. The same line 


* Private communication. 

7 By a method similar to that described by A. G. Nickle, 
Dissertation, Harvard University (1941). 

* H. Gershinowitz and E. B. Wilson, Jr., J. Chem. Phys. 
6, 247 (1938). 

** A search of the literature failed to reveal any infra- 
red data on this compound, but positions of Raman lines 
hs Kahovec and K. W. F. Kohlrausch, Zeits. f. physik. 
Chemie 46, 165 (1940) ] do not correspond to the observed 
absorption maxima. 

t Since no infra-red data was found in the literature and 
the Raman spectrum [Lespieau, Bourguel, and Wakeman, 
Bull. Soc. Chim. Fr. 51, 400 (1939) ] indicated the possibil- 
ity that lines observed in our material might be due to 
methyl cyclopropane, some was prepared [N. Demjanoff, 
21 (1895)] and a tentative absorption spectrum 

ained. 


1928)" Coffin and O. Maass, J. Am. Chem. Soc. 50, 1427 
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may appear in the Raman effect since that was 
studied using the sample first prepared, after the 
preliminary investigation of the infra-red. 


Spectra 


The automatic prism spectrometer used in this 
series of investigations in the infra-red has been 
described previously.!° The instrument was cali- 
brated by comparison of known bands of COs, 
CsHy, and CHs3Br with the calculated 
calibration curve of the prisms used. The CaF» 
prism was used in the range 2 to 8.54, the NaCl 
prism in the range 8 to 144, and the KBr prism 
in the range 13 to 25y. 

The cells used were the 30-cm Pyrex tubes, 
with KBr windows sealed on the ends with clear 
Glyptal, which have been used in previous in- 
vestigations. To prevent attack on the stopcock 
grease by the hydrocarbon a stopcock with 
capillary leads was used, and a bead of Hg held in 
the capillary by the pressure inside the cell kept 
the gas from coming in contact with the grease. 

Figure 1 gives the observed spectrum of 
cyclobutane. Results of runs in which the sample 
was contaminated with butene-1 are not included. 
Everywhere, except for the range above 3200 
cm~ the spectrum has been covered more than 
once, and all bands reported in Table I except 
those marked with ? were observed in all runs at 
the higher pressures. 

In the spectroscopic investigation of impurities 
(see above) special attention was paid to the 
bands at 754 and 737 cm~. These showed no 
change in intensity within the limits of error of 
measurement (about 2 percent on the percent 
transmission scale) after a distillation which 
changed the intensity of a comparable line of 
butene-1 by a factor of 2. Bromination had no 
effect on their intensity either. If they are due to 
an impurity its properties must, then, be very 
similar to those of cyclobutane and we have been 
unable to find one satisfying the physical require- 
ments which cannot be eliminated on the basis of 
spectroscopic evidence, or a consideration of the 
starting material. The bands mentioned have, 
therefore, been assumed to be actually due to 


10 H. Gershinowitz and E. B. Wilson, Jr., J. Chem. Phys. 
6, 197 (1938) and A. J. Wells, Dissertation, Harvard Uni- 
versity (1941). Latter describes changes made since pub- 
lication of I, including addition of CaF: prism. 
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PERCENT TRANSMISSION 


7800 2000 7400 2600 7800 3000 3200 3700 7200 4700 


FREQUENCY IN WAVE NUMBERS 


Fic. 1. The infra-red absorption spectrum of cyclobutane. Pressures: (a) 165 mm; 
(6) 135 mm; (c) 42 mm; (d) 3 mm; (e) ca. 1 mm of Hg. 


cyclobutane, which necessitates postulation of an 
unobserved branch, probably on the low fre- 
quency side of the observed pair. 

The Raman spectrograph used has also been 
described" previously. Qualitative polarization 
measurements were made by the “Polaroid” 
method.” The use of dyes and filters to eliminate 
the Hg green lines and excitation from the 4047 
Hg line, which has been described in other papers 
in this series,“ proved especially helpful in 
dealing with this material. 

The Raman spectrum of liquid cyclobutane at 
25°C was obtained. The small quantity of cyclo- 
butane prepared necessitated the use of a par- 
ticularly small Raman tube (capacity20.5 cc) 
which is described elsewhere." As shown by the 
microphotometer traces (Fig. 2) the spectra ob- 
tained, while satisfactory in most respects, might 
not show up very weak shifts since the intensity 
of the continuous background is rather higher 
than usual in comparison with that of the light 
from Raman scattering. The traces also show the 
difficulty in separating lines. Careful examination 


11 H. Gershinowitz and E. B. Wilson, Jr., J. Chem. Phys. 
6, 247 (1938). 

2 J. T. Edsall and E. B. Wilson, Jr., J. Chem. Phys. 6, 
124 (1938). 

13Cf, A. J. Wells and E. B. Wilson, Jr., J. Chem. Phys. 
9, 314 (1941). 

4H. D. Noether, Dissertation, Harvard University 
(1942). 


of plates and microphotometer traces has led us 
to the conclusion that the shifts centered around 
_ 1450, 1220, and 933 cm are probably all 
doublets, and that the broad band around 3000 
cm! is composed of three lines, but the spacing 
is so small that exact determination of the posi- 
tions of the components is impossible. All shifts 
reported may be out by at least +5 cm™. 

The very weak line reported at 750 cm~ was 
observed on two plates though its intensity was 
so low that its separation from the background 
was not definite on the microphotometer traces. 
It cannot be regarded, then, as any more than 
probable, and if it is present we believe it to be 
due to an impurity. This would correspond to the 
728 cm— band observed in the infra-red spectrum 
of the sample of cyclobutane used for the Raman 
work (see section on purity of product). 

Table I gives the frequencies of the absorption 
maxima observed in the infra-red spectrum and 
their relative intensities, together with the proba- 
ble band type (perpendicular or parallel) and an 
assignment of ‘‘P,” “Q,” and ‘‘R” branches in 
cases where these could be determined (see 
assignment). The table also includes shifts ob- 
served in the Raman effect, the exciting lines 
from which these shifts were observed, relative 
intensities, and an indication of whether or not 
they appear polarized. 
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THEORETICAL 
Symmetry of Cyclobutane 


Recent calculations by Aston et al'® on the 
configuration of cyclopentane have presented the 
possibility that the ring in that substance is 
deformed in such a way that one of the carbon 
. atoms is out of the plane of the other four. If so, 
then one might be justified in considering a non- 
planar structure of the ring in cyclobutane as 
well. 

The symmetry of cyclobutane if the ring is 
distorted out of a plane is that of the point group 
Co», whereas if the ring is planar the symmetry 
is Dy. If it is the latter a center of symmetry is 
present and in such a case the selection rules 
forbid the appearance of frequencies corre- 
sponding to the same vibration in both the 
Raman and infra-red spectra. If, on the other 
hand, the symmetry is C2,, all vibrations active 
in one may also appear in the other. In the ob- 
served spectra some approximate coincidences do 
occur (at 3000, 1450, 1240 and possibly 900 cm-') 
but all these can be explained as readily on the 
basis of the D4, as on the C2, selection rules. 
Accidental coincidences due to similarity, but not 
identity, of types of molecular vibrations are ex- 
pected, and bands at the frequencies mentioned 
occur in most hydrocarbons. The points of greater 
interest are the cases where lines occurring in 
one spectrum do not appear in the other. The 
band in the infra-red spectrum at 627 cm has 
apparently no analogue in the Raman spectrum 
and the line at 1010 cm™ in the latter almost 
certainly does not appear in the former. Since the 
latter line is highly polarized it is probably the 
symmetrical ‘‘breathing’” frequency of the ring 
(see assignment) and this should appear in the 
infra-red spectrum if not forbidden by selection 
rules. 

Besides the lack of coincidence in the two 
spectra we may point out the small number of 
bands appearing in the infra-red. Only four strong 
bands are observed in the range 400-1500 cm-', 
in agreement with predictions from the Da 
selection rules, but definitely not in agreement 
with those of C2, which would permit the appear- 


6 J.G. Aston, S. C. Schumann, H. L. Fink, and P. M. 
Doty, J. Am. Chem. Soc. 63, 2029 (1941). 


373 


ance of sixteen bands. The same argument dis- 
poses of the possibility that the symmetry is V, 
(distortion of the ring in-the plane), in which case 
eight fundamentals should appear. Since funda- 
mentals permitted by the selection rules to 
appear in the infra-red usually do so, the con- 
clusion seems obvious. 

Table II gives selection rules and polarization 


TABLE I. 


Infra-red 


Intensity 
and¢ 
polarization 


Intensity 
Frequency and 
(cm~!) types 


4500 
4250 


2988 
2962 (e, k) 
2938 


2880 (e, k) 


0 
8 
4 
1 
2 
3 
0 
4 
5 
0 
1 
0 
2 
0 
1 


~1110? (e)? 
1010 (e, f, 


2 

3 (?) 
0 
10 (P) 


933 (9995) 7b 
~ 750? (e)? 0 (P?) 


R 
4 5 (1) 


@ See text for differentiation between types of bands. 

* The letters in parentheses following the frequencies of the nea 
shifts indicate the exciting lines from which these shifts were o 
Kohlrausch’s nomenclature is employed. 


Raman 
3950 
3775 
3350 
3190 
8 
2990 2 12 
2900 8 (P) 
2735 
ind 2626 3 (P) 
2463 
all 3360 
1 
ing 2145 
2090? 
2050 
ifts 1987? 
1896 
1827 
vas 1653 Rt 
1636 Q? 
was 1607? P? 0 
ind 1565 R? 
1547. Q? 1 
ces. 1528 P? 
lan 1511? 0 
1466 R 1465 4 (P) 
be 1453 Q 10(1) | 1445 ©) 4 
the 1436 
1371 R 
um 1355 2 (1?) 
nan ~1340 
1987 10(.L) 
1 4: 
1243 
10n 1216 
and 1 
1070? 0 
)ba- 925 R 
902 ( 
in 
in 
ob- 644 | 
ines 
tive 
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or band types of the active vibrations of cyclo- 
butane if the symmetry is D4,.15¢ 


Assignment of Frequencies 


A description of the types of vibrations oc- 
curring in the various representations is given in 
Table III. We are using the nomenclature de- 
veloped by Linnett'® to differentiate between the 
various possible motions of the methylene group, 
except that the motion he referred to a bending 


4358 A 
Hg group 


Fic. 2. Microphotometer traces of Raman spectra of 
cyclobutane. All plates taken with 0.09 g/l. Rhodamine 5 
GDN extra in the alcohol circulated in the cooling system. 
Filters used besides this were: (A) Wratten No. 2A around 
the Raman tube; (B) no other filters; (C) Polaroid around 
the Raman tube to polarize incident light parallel to the 
axis of the tube. 


we prefer to call wagging both to prevent con- 
fusion with the so-called bending motions of the 
CH; group which are more closely related to the 
“CH deformation’ and because we feel the 
latter term to be somewhat more descriptive. 
This system we have found quite convenient 
both in discussions and in calculations. 

Starting with the simplest and perhaps most 
definite case, we may assign frequencies to the A’ 
vibrations on the basis of the Raman polarization. 
It will be seen from Table I that only four lines 
show definite polarization. The 2940 and 2880 


15a J. E. Rosenthal and G. M. Murphy, Rev. Mod. Phys. 
8, 317 (1936). 
16 J. W. Linnett, J. Chem. Phys. 6, 695 (1938). 
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cm lines are clearly the symmetrical C—H 
stretching and 21450 resonating with it. 1465 
cm is the CH: deformation and 1010 the 
symmetrical ‘‘breathing”’ frequency of the ring 
by obvious frequency arguments. The latter is 
close to a prediction by Aubert!’ but does not 
check too well with what Kohlrausch and 
Skrabal'® expected from work on derivatives. 

The other representations on which we have 
sufficient evidence to make a definite assignment 
are the infra-red active A” and E”. In cyclo- 
butane it is easily possible to differentiate be- 
tween vibrations with these symmetries by the 
band envelopes. Calculations from the moment 
of inertia of the molecule using Gerhard and 
Dennison’s'® formula for the P—R spacing of 
parallel bands in symmetrical top molecules lead 
us to expect this spacing to be 39 cm for A.” 
bands. We observe spacings of 40 cm in the 
902 cm band, and 30, 29, and 27 cm, re- 
spectively, in the 1453, 1257, and 627 bands. 
Eliminating 737 cm as a possible fundamental 
(this will be considered later), it is obvious that 
902 cm- is the A2” rocking* and the other three 
are the E”’ fundamentals. The normal coordinate 
treatment in E” shows that 1453 cm~ is princi- 
pally CH, deformation, 1257 cm CHe wagging, 
and 627 cm the degenerate ring deformation. 
The low frequency of the ring deformation we 
believe to be primarily due to the large inter- 
action between the C—C stretching and the 
C—C—C angle bending (see Table V). The 
C—H stretchings are uncertain but surely near 
2990 cm~ in both A2”’ and E”’. The values given 
in Table III for these motions are guesses. 

The other directly observed frequencies are the 
depolarized Raman shifts. Of these, that at 1445 
must be the Be’ CH: deformation since it is 
the only line left near the expected frequency. In 
accordance with Kohlrausch’s!* predictions, and 
because the line is the strongest in the proper 
frequency range, we assign the 940 cm shift to 
the unsymmetrical ring stretching in B,’. The 
unpolarized CH stretches at 2988 and 2962 are 

17M. Aubert, Comptes rendus 203, 661 (1936). 

18 K, W. F. Kohlrausch and R. Skrabal, Zeits. f. Elec- 
trochemie 43, 282 (1937); or K. W. F. Kohlrausch, Der 
Smekal-Raman effect, Erganzungsband (Springer, Berlin). 
(1933) L. Gerhard and D. M. Dennison, Phys. Rev. 43, 197 


*See also Table III for an a priori calculation of A.” 
frequency. 
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assigned as indicated in Table III because the 
form of the force constants in the Raman active 
representations leads us to expect the E’ fre- 
quency to have a value between those in Aj’ 
and B,’. 

The remaining lines observed are not so readily 
or definitely assignable. To distinguish between 
them we have resorted to the results of the 
normal coordinate treatment. Using Stitt’s 
value”® for the C—C—H bending constant, Hz, 
and a value of the C—C—C bending constant 
estimated from Wilson and Wells’ work” on 
propylene we have come to the conclusion that 
the frequencies observed at 1228 cm, 1209 cm“, 
and 923 cm are probably to be assigned as 
indicated in Table III. These calculations also 
lead to a prediction of about 600 cm“ for the B,’ 
in-plane ring bending. A careful examination of 
this region on the Raman plates has failed to 
reveal any line below 750 cm. Since an alterna- 
tive explanation of the 750 cm“ line is available 
(see under Spectra) and since it is definitely 
higher than the value we expect for the funda- 


TABLE II. 


Symmetry with respect to 
Number of 
funda- 
mentals 


Representa- Selectiont 


tion 


~ 
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rules, but such a violation might occur in the 
liquid phase spectrum. The appearance of fairly 
strong lines at very nearly this frequency in 
all cyclobutane derivatives investigated by 
Kohlrausch, Reitz, and Skrabal!* ” (except cyclo- 
butanone where the strain introduced by the 
carbonyl group may account for its absence) 
serves to confirm our hypothesis. We have 
assigned it to A,’ on the basis of calculations 
from the constants used in E’. 


TABLE III. 


Frequency 


Repre- 


sentation Type of vibration Assigned* 


C —-H stretching 
deformation 
ring stretching . 
wagging 
CH2 wagging 
ring stretching 

C —H stretching 
CHe deformation 
ring bending (in plane) 
C —H stretching 
CH: twisting 
CHe rocking 
CH: twisting 

C —H stretching 
CHe rocking 

C —H stretching 
CH: rocking 

ring bending (out of plane) 
CHe twisting 

C —H stretching 
CH2 deformation 
CHe wagging 
ring deformation 


2940, 2880 
1465 


A 
degenerate 


| 


degenerate 


* Planes | plane of the ring and bisecting the C —C —C angles. 
R=Raman active, J—R=infra-red active, F =inactive in both 
spectra, P =polarized, dp =depolarized, || =parallel type band, | =per- 
pendicular type band. 


mental, we are of the opinion that the shift 
corresponding to the fundamental was not 
observed. 

The doubtful Raman line at 1110 cm~, if 
present, we believe to be one of the non-de- 
generate CH, twisting frequencies. The appear- 
ance of either would violate the D«, selection 


20F. Stitt, J. Chem. Phys. 7, 297 (1939). 
*1E. B. Wilson, Jr., and A. J. Wells, J. Chem. Phys. 9, 
319 (1941). 


See text. 


The fundamentals not observed we have as- 
signed from consideration of the combination 
bands appearing in the infra-red spectrum, to- 
gether with calculations using Stitt’s H,?° and 
consideration of the spectrum of derivatives.'* ” 
The two ring bending motions, in and out of the 
plane, were assigned on the basis of the two 
lowest definite combination bands appearing in 
the infra-red. The 737 cm band is believed to be 
a combination since it may reasonably be ex- 
plained as such as indicated in Table IV, whereas 
627 cm= could not be. Also, its intensity is 
probably low for a fundamental but high for a 
difference ‘‘tone,’’ the other possible source of 
such a band if impurities are excluded. Unfortu- 
nately this band is so weak that we have been 


2K. W. F. Kohlrausch and R. Skrabal, Monats. f. 
Chem. 70, 44 (1937); O. Reitz and R. Skrabal, Monats. 
f. Chem. 70, 398 (1937); R. Skrabal, Monats. f. Chem. 70, 
460 (1937). 
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unable to determine its envelope, which would 
confirm or disprove the assignment in Table IV. 

It will be noted that the assignment of funda- 
mentals given in Table III offers an explanation 
of all bands observed with the exception of the 
doubtful one at 1511 cm~. That is very weak 
and in the region of strong absorption by 
atmospheric water so its presence is distinctly 
uncertain. 


Calculations 


In all calculations on cyclobutane we have 
assumed, in the absence of experimental evidence, 


TABLE IV.* 


Observed 


frequency Assignment 


737 

~1070? 

1216 

1355 

1511? 

1547 

1565 

1607? 

1636 


1827 


1896 


1987? 
2050 
~2090? 
2145 
2183 
2272 


145+ 595= 740° 
145+ 923=1068 
595+ 627 =1222 
145+1209 = 


627+ 923=1550 
627+ 940=1567 
145+ 1445 =1590 
627+ 1010 = 1637 
902+ 923=1825 
627+ 1209 = 1836 
902+ 1010 =1912 
627+ 1260 = 1887 
940+ 1050 = 1990 
595+ 1453 = 2048 
627+ 1465 = 2092 


2900+ 1200 


1050+ 1230 =2280 
1110+1260 = 2370 
902+ 1465 = 2367 
923+ 1453 =2376 
1010+ 1453 = 2463 
1209+-1257 = 2466 
1200+ 1445 = 2645 
1257+1465 =2722 
1260+1453 =2713 
2x 1450 =2900 
1453+1465 = 
1445+ 1453 =2898 
1010+2900+1200 (seeabove) A,.’’, E” 
1010+ 2360 (see above) 
3X 1257 =3771 
other combinations of three 
v’s 
1010+1257+1465 = 3732, etc. 
1010+2980 = 3990 
940+ 2990 = 3930, etc. 
2~3000+1200 Ag’ 
3000+ 1450 


* This table certainly does not include all combinations which might 
explain the observed lines. Except in a few cases where it was neces- 
sary, combinations of three fundamentals are not used to explain 
lines. These undoubtedly could be observed in a case such as this 
where the frequencies occur in groups. 


@ See text. 
(R) Observed in the Raman spectrum. 
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that it has the following dimensions: 


C—C=1.53A,* 
C—H=1.09A,™ 
C—C-—C angle=90°, 
H—C-—H angle=109° 28’ (tetrahedral), 
C—C-—H angle=114° 6’. 


Another plausible assumption would have been 
to take the C—C—H and H—C—H angles equal 
(112° 59’). The difference is negligible for our 
purposes. With the same force constants, fre- 
quencies are different by less than 2 percent 
when calculated on the two assumptions. 

The secular equation was set up in valence type 
symmetry coordinates*®?* and the results were 
checked by repeating the process using a different 
method recently developed by C. W. Lippman.?? 
We would like to take this opportunity to thank 
the latter for giving us the details of his method 
prior to its publication. C—H stretches were 
factored out”® and treated separately. The treat- 
ment showed Stitt’s Ky? to be a satisfactory 
force constant for this molecule. 

The internal coordinates selected for treatment 
of this problem were the sets of changes in the 
four C—C distances, the eight C—H distances, 
the four H—C—H angles, the sixteen C—C—H 
angles and the four C—C—C angles. The sym- 
bols used here to represent these changes are, 
respectively, R;, 7:;, a:, B:, and y;. In Fig. 3 are 
shown some of the coordinates around a single 
carbon atom to illustrate the system of number- 
ing coordinates. Subscripts on a’s and y’s, not 
shown, are the numbers of the central carbon 
atom of the corresponding angles. 

The out of plane bending of the ring was not 
treated and another type of coordinate would 
have to be added to include it. Omission of this 
motion does not necessarily imply that the By,” 
representation, in which it occurs, cannot be 
treated since it is such a low frequency that 


23. Pauling and L. O. Brockway, J. Am. Chem. Soc. 59, 
1223 (1937). 

* N. Ginsburg and E. F. Barker, J. Chem. Phys. 3, 668 
(1935). 

25 E, B. Wilson, Jr., J. Chem. Phys. (a) 7, 1047 (1939) ; 
(b) 9, 76 (1941). 

.%M. A. El’yashevich, J. Phys. Chem. U.S.S.R. (a) 14, 
1381 (1940); (b) 15, 831 (1941). 

27 C, W. Lippman, soon to be published. 


| 
Sym- 
metry 
A 
E” 
E” 
A 
E” 
A 
A 
_ 
E" 
2360 
2463 
2626 
2735 
2880 (R) 
2900 
3190 
3350 
3775 
3950 
4250 
4500 
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interactions with other vibrations in the repre- 
sentation are probably negligible. If so, it may 
be factored out by setting the corresponding 
force constant equal to zero. 

The four types of motions of the CH» group 
(see above) may be represented in terms of 6’s 
(for carbon atom C;) as: 


CH: deformation = 
CHe wagging = 
CH, twisting 

rocking = 8; 


Analogous combinations represent the same mo- 
tions on other carbon atoms. 

The redundancy conditions on the coordinates 
used are: 


(1) At each carbon atom: 

2V/ 2a: + +27: =0; 
(2) 
(3) 2(Ri—Rs) + 4) =0; 
(4) 2(R2—Rs) — +4) =0. 


A complete treatment of the vibrational prob- 
lem has proved possible assuming Stitt’s?® con- 
stant Hg for the C—C—H angle change. As 
shown in Table III the maximum deviation, 
+3.2 percent of calculated from definitely ob- 
served frequencies, occurs in E’’, where the 
number of interactions possible but neglected 
is larger than in the other cases. The fit could be 
improved by adjusting some of these, but in view 
of the uncertainties both in the actual value of 
the C—C—H angle constant and in the geome- 
try of the model, it was felt that this would only 
add to the complexity of the potential function 
without necessarily improving it. The C-C—H 
angle constant has not been evaluated from the 
data observed here since a number of interactions 
between these coordinates appear to be of some 
importance. On the average, Stitt’s value appears 
to be as satisfactory as possible without including 
these interactions. 

The force constants used in the potential 
function are given in Table V. The constants 
given are for internal coordinates except for Fa 


*The “CHe deformation” as described b 
reference 16, requires a change in the a as wal. 
on the carbon atom. Actuall 
corresponding + are involved in a redundancy relation and 
sO are not independent. 


Linnett, 
as the 
these coordinates and the 
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and F(Z’). The former is the constant for the 
combination of changes in the angles around one 
carbon giving rise to the ““CH» deformation’”’ 
when the corresponding element in the G matrix”® 
is obtained from a only. Due to the redundancy, 


Fic. 3. Here the actual bonds and angles are labeled 
R, r, etc. In the text it is the corresponding changes to 
which these symbols refer. 


fa always appears in this combination. If we 
assume Stitt’s H,?° to be transferable, Fz may be 
broken down, leading to a value of fag=0.047. 
This is approximately four times the value Stitt 
found in ethane, but still possible. If, however, we 
assume fag=0 then f.=0.428 (cf. 0.547 from 
Stitt’s H.) which is an unexpectedly large 
change. The truth may well lie somewhere be- 
tween these two. 

Fgg(E’) is the interaction appearing in the E’ 
representation between the two combinations of 
6’s appearing there. This is the only interaction 
between angles which we can be nearly certain is 
present, although in view of the uncertainty in fg, 
its numerical value is not definite. 


Discussion of Force Constants 


The only force constant evaluated which shows 
definitely a change from that obtained in 
unstrained cases is fr. This is reduced about 5 
percent from the value obtained by Stitt®® in his 
work on ethane, and about 8 percent from the 
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TABLE V. 


Coordinates 


Value involved 


dynes 
c 


Sr 4.26X 105 Ri 


dyne 
radian 

dyne 
radian 

dyne 
radian 


—0.059 x 10-3 the type Riai 


f Ra 
Srp 


Sry —fry 


0.108 x 10-8 the type Rif; 


1.19 «10-3 the types 
Riyst 
dyne-cm 
radian? 
dyne-cm 
radian? 
dyne-cm 
radian? 
dynes 
cm’ radian? 


F,** 0.690 


fp (reference 20) 0.654 B; 


Fop(E’) — 0.023 x 


* The interaction constants given here are taken from matrix expres- 
sions for the potential energy. For use in the quadratic form of this 
quantity they none be multiplied by 2. 

** Fa =fq —4(2/5) +(2/5)fp. 
text. 

In that order. . 

In both degenerate representations there are, of course, two equiva- 
lent, orthogonal sets of coordinates, of which only one was considered 
in the calculations. In each set two combinations of #’s appear as 
symmetry coordinates. These have been denoted by the superscripts ’ 
and ”’, In E’ the two 8 coordinates in the set used are 


B'(E’) =}Roro(ti +t2—ts —ts), 

B’'(E’) =} Roro(v1 —v2 —03 +4), 
in the notation given in the text. 
value predicted for a C—C bond length of 1.53A 
by Crawford and Brinkley’s*’ relation. This is in 


28 B. L. Crawford and S. R. Brinkley, Jr., J. Chem. Phys. 
9, 69 (1941). 


qualitative agreement with results of Kohlrausch 
and co-workers.” 

We have found no suitable constant with which 
to compare f, so determination of the effect of 
strain on this will have to wait until one is 
evaluated. Since the corresponding angle is con- 
siderably strained we might expect an appreciable 
effect. 

Two other interesting features of the potential 
function obtained, presumably effects of the 
strained configuration, are the absence of inter- 
actions between C—C stretches and the large 
interaction between R and y. Also, the agreement 
between Fra* and Stitt’s Ke,, analogous con- 
stants corresponding to the interaction between 
R and a symmetrical deformation of the angles 
around the carbons which it connects, is rather 
striking. Unfortunately, until further evidence is 
available on hydrocarbons, both cyclic and non- 
cyclic, it appears useless to try to develop general 
theories on their potential functions. 

I should like to take this opportunity to ex- 
press my thanks to Professor E. B. Wilson, Jr., 
both for suggesting this problem to me and for 
his continued interest and many helpful sug- 
gestions throughout the course of the work. 


29 Kohlrausch and Seka, Ber. 69, 729 (1936); Kohlrausch 
and Skrabal, Monats. f. Chem. 70, 44 (1937). 

*Interaction between the CH: deformation and the 
C—C stretching symmetry coordinates. 


THE JOURNAL OF CHEMICAL PHYSICS 


VOLUME 11, NUMBER 8 AUGUST, 1943 


Some Remarks on the Kinetics of Vinyl Derivative Polymerization 


J. Apere, G. Go_princer, H. MARK, AND H. Narpus 
Queens College, Flushing, New York and Brooklyn Polytechnic Institute, Brooklyn, New York 


(Received March 4, 1943) 


The various elementary steps of vinyl polymerization 
reactions are discussed. Activation is considered to occur by 
collision of an unactivated monomer either with another 
monomer, a catalyst, a solvent, or a polymer molecule. It is 
further taken into account that certain types of active 
nuclei can be more easily deactivated than others. Two 
steps of propagation are considered: normal chain growth 
and chain transfer. Termination of growing chains is 
assumed to be possible either by collision with another 
growing chain, with an activated monomer, with a mole- 
cule of the catalyst, or with certain chain breaking im- 
purities. Using the assumption that for a certain period of 


the reaction a steady state of activated centers is main- 
tained, equations are derived for the initial over-all rate of 
the reaction and for the number average polymerization 
degree. These equations are compared with experiments 
which have been carried out on the polymerization of 
styrene at various temperatures, monomer concentrations, 
and catalyst concentrations. Comparison between theory 
and experiments shows that in first approximation, the 
initial over-all rate of the reaction ana the initial number 
average polymerization degree can be represented by the 
derived equations. 


I. INTRODUCTION 


N the last few years a number of experimental 

studies on the kinetics of vinyl-type polymer- 
ization have been published, investigating the 
rate of these processes as functions of various 
experimental parameters, such as chemical na- 
ture, concentration (or better activity) or partial 
pressure of monomers, temperature, chemical 
character and amount of catalyst, nature of 
solvents, etc.! 

The monomers used were mainly styrene, 
vinyl acetate, vinyl chloride, and acrylic acid 
derivatives, the catalysts peroxides, ozonides, 
free radicals or metal halides, the solvents 
hydrocarbons, ketones, substituted phenols, and 
chlorinated hydrocarbons. The temperature 
range was between 60° and 150°C.? 


1R. E. Burk, H. E. Thompson, A. J. Weith, and I. 
Williams, Polymerization (Reinhold Publishing Company, 
New York, 1937); H. Mark and R. Raff, High Polymeric 
Reactions (Interscience Publishers, New York, 1941). 

2J. Abere, G. Goldfinger, H. Mark, and H. Naidus, 
Paper presented at the Buffalo Meeting of the American 
Chemical Society in September, 1942; H. N. Alyea, J. J. 
Gartland, and H. R. Graham, Ind. Eng. Chem. 34, 458 
(1942); H. Dostal, H. Mark, and R. Raff, Ind. Eng. Chem. 
29, 595 (1937); A. C. Cuthberson, G. Gee, and E. K. 
Rideal, Nature 140, 889 (1937); C. S. Marvel and J. C. 
Cowan, J. Am. Chem. Soc. 61, 3156 (1939); C. S. Marvel, 
J. Dee, and H. G. Cooke, J. Am. Chem. Soc. 62, 3499 
(1940); C. S. Marvel, G. D. Jones, T. W. Mastin, and G. L. 
Schertz, J. Am. Chem. Soc. 64, 2356 (1942); H. W. Mel- 
ville, Proc. Roy. Soc. A165, 163, 511 (1937); A167, 99 
(1938); J. K. Moore, R. E. Burk, and H. P. Lankelma, J. 
Am. Chem. Soc. 63, 2954 (1941); R. G. W. Norrish and 
E. F. Brockman, Proc. Roy. Soc. A163, 205, (1937); A171, 
147 (1939) ; Trans. Faraday Soc. 35, 1087 (1939); R. G. W. 
Norrish and R. R. Smith, Nature 150, 336 (1942); H. 
Pfann, Thesis, Brooklyn Polytechnic Institute (May, 


The results of these investigations provide a 
very useful background for developing a quanti- 
tative scheme for such composite processes. 
All investigators agree that they involve a 
number of characteristic individual steps, such 
as activation and deactivation of the monomer, 
propagation of the activated monomer into long 
chains with an activated end (or two activated 
ends), transfer of the activation to another 
molecule (monomer or solvent), branching of the 
growing chains and termination of the growth in 
various ways.’ The first task for a numerical 
treatment of such composite reactions is always 
to find out how the different elementary steps 
cooperate in producing and consuming various 
intermediate (short-living) configurations and 
how they finally build up the stable end-product 
with its characteristic properties, such as aver- 
age polymerization degree, molecular size dis- 
tribution, frequency of branching, etc. In prac- 
tically all of the papers quoted above, the authors 
have carried out an analysis of this kind and in 
most cases succeeded in deriving formulas which 
closely agreed with their experimental findings. 
1942); C. C. Price and R. W. Kell, J. Am. Chem. Soc. 63, 
2798 (1941); C. C. Price, R. W. Kell, and E. Krebs, J. Am. 
Chem. Soc. 64, 1103 (1942); C. C. Price and D. A. Durham, 
J. Am. Chem. Soc. 64, 2508 (1942); G. V. Schulz, H. 
Dinglinger, and E. Husemann, Zeits. f. physik. Chemie 
B34, 187 (1936); B36, 184 (1937); B39, 246 (1938); B43, 
385 (1939) ; B50, 305 (1941) ; Naturwiss. 27, 387, 659 (1939) ; 
H. Staudinger and G. V. Schulz, Ber. 68, 2320 (1935). 

’P, J. Flory, J. Am. Chem. Soc. 59, 241 (1937); R. 
Ginell and R. Simha, J. Am. Chem. Soc. 65, 706, 715 
(1943); H. Mark and R. Raff, High Polymeric Reactions 


(Interscience Publishers, Inc., New York, 1941); New York 
Acad. Sci., Discussion on High Polymers (January, 1943). 
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As necessary, however, in pioneering work of 
this kind, emphasis was mainly laid on the special 
case under consideration and the experimental 
parameters listed above were not always varied 
over as wide a range as would be desirable for a 
systematic study of these processes. It was there- 
fore felt that it might be appropriate to help 
extend our present knowledge by a more sys- 
tematic study of a few particularly interesting 
systems. Two monomers (styrene and metha- 
crylic ester) were selected and their polymeriza- 
tion was (and is still being) studied in several 
typical solvents (methyl ethyl ketone, carbon 
tetrachloride, toluene, benzene) over a fairly 
wide range of catalyst concentration using 
different types of catalysts (benzoylperoxide, 
p-Br-benzoylperoxide, tin tetrachloride), various 
initial monomer concentrations (10-80 percent), 
various types and activities of inhibitors, and 
working at temperatures between 60 and 150°C. 
Some of the results have been presented previ- 
ously and are being published in more detail 
elsewhere.‘ In order to evaluate and coordinate 
them it seemed useful also to extend somewhat 
the formal kinetic scheme for such processes. This 
is the purpose of the present paper. 

It must be emphasized, however, that formal 
kinetics is not the final goal of the quantitative 
study of a polymerization reaction. What we 
would like to arrive at is a full understanding of 
the chemical nature of each single elementary step 
such as activation, propagation, termination, 
etc. To reach this goal, however, it seems neces- 
sary first to describe each of these steps quantita- 
tively and, if possible, to assign to each of them 
a certain order of reaction, a definite frequency 
factor and a specific activation energy. Such 
numerical characterization is very helpful in 
narrowing down the speculations as to the chem- 
ical nature of these individual elementary proc- 
esses, and eventually allows the interpretation 
of them as a certain change in the electronic 
structure of the molecules involved (opening and 
closing of covalent bonds, polarization of the 
double bond, displacement of a pair of elec- 
trons, etc.). 

4J. Abere, G. Goldfinger, H. Mark, and H. Naidus, 
Paper presented at the Buffalo Meeting of the American 


Chemical Society in September, 1942; New York Academy 
_— Discussion on High Polymers, January 8 and 9, 


Systematic formal treatments of polyreactions 
have been carried out previously by Breiten- 
bach,® Flory,® Raff,’ and Schulz’ and the general 
aspect of vinyl polymerization has been qualita- 
tively correlated with the fundamental proper- 
ties of the aliphatic double bond. This correla- 
tion has been very definitely improved and ex- 
tended to other characteristic properties of the 
double bond (cis-trans-isomerization, addition 
reaction) by the recent work of Eyring and his 
collaborators.°® 

In all the following calculations, we shall use 
the Bodenstein principle, assuming that after a 
comparatively short time a steady state of 
active (growing) centers is established and main- 
tained throughout a considerable period of the 
experiment. This simplifying assumption has 
been always used, and is presumably in most 
cases (moderate amount of polymer being 
formed, presence of solvents) good enough for 
the interpretation of the experiments. It must, 
however, be kept in mind that it represents an 
oversimplification of the actual conditions, and 
it is of importance that Ginell and Simha* have 
recently carried through the integration of the 
general differential equations of a polymeriza- 
tion reaction without any restricting assump- 
tion of this kind. Their results are of great value 
to decide whether under given conditions the 
Bodenstein principle can be legitimately applied 
or not. 


II. DISCUSSION OF THE ELEMENTARY REACTIONS 


Let us first consider the activation or nucleus 
formation. It is not yet definitely established 
whether pure thermal activation can take place 


5]. W. Breitenbach and R. Raff, Monats. f. Chem. 69, 
1107 (1936); J. W. Breitenbach and H. Rudorfer, Monats. 
f. Chem. 70, 37 (1937); J. W. Breitenbach, Monats. f. 
Chem. 71, 275 (1938); J. W. Breitenbach and H. Maschin, 
Zeits. f. physik. Chemie A187, 175 (1940). 

6P. J. Flory, J. Am. Chem. Soc. 59, 241 (1937); 62, 
1057, 1561, 2255 (1940). 

7H. Dostal, H. Mark, and R. Raff, Ind. Eng. Chem. 29, 
595 (1937); H. Mark and R. Raff, High Polymeric Reactions 
(Interscience Publishers, New York, 1941). 

8G. V. Schulz, H. Dinglinger, and E. Husemann, Zeits. 
f. physik. Chemie B34, 187 (1936); B36, 184 (1937); B39, 
246 (1938); B43, 385 (1939); B50, 305 (1941); Naturwiss. 
27, 387, 659 (1939); H. Staudinger and G. V. Schulz, Ber. 
68, 2320 (1935). 

®R. A. Harman and H. Eyring, J. Chem. Phys. 10, 557 
(1942); J. L. Magee, W. Shand, and H. Eyring, J. Am. 
Chem. Soc. 63, 677 (1941); R. G. W. Norrish and R. R. 
Smith, Nature 150, 336 (1942); T. Ri and H. Eyring, 
J. Chem. Phys. 8, 433 (1940). 
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by the collision of two monomers, but the course 
of the polymerization of ethylene and styrene 
at high temperatures in the gaseous phase seems 
to indicate its existence. The possible results of 
such a step are :* 


CH,.=CH2+CH2= 
CH.=CH:+ —CH.—CH:2-, (1) 
—CH,—CH.—CH2—CH:2-, (2) 
(3) (1) 
—CH,—CH;+—CH=CH,, (4) 
CH.=CH.+CH:=CH2— 
CH;—CH;+ —CH=CH-. (5) 


Cases (1u,1) to (1u,3) represent only ex- 
change of energy between the colliding particles, 
while (iu,4) involves the exchange of one 
hydrogen atom and (1u,5) the exchange of 
even two of them. The activation energy for 
these two cases is presumably rather high'® and 
they will be infrequent at temperatures below 
200°C. The third reaction involves also a high 
activation energy and—even if it occurred—will 
presumably have the tendency to go over into 
(1u, 2). Hence, we shall only include (1u, 1) and 
(1u, 2) in the further discussion. In the following 
we shall always denote the monomer with my, 
the dimer with mo, etc.; activation on one or 
both ends of the chain is indicated by an asterisk. 
The two above cases can therefore be written as 


(1u, 1) 
(1, 2) 


Their nature as activation processes is indi- 
cated by the fact that there are no asterisks on 
the left-hand sides, but there are some on the 
right-hand sides of the equations; activation is 
produced. 

Another process of nucleus formation is the 
collision of m; with a catalyst c, which take 


m,+m,—*m2*. 


_*We shall use the following symbols to distinguish 
different reaction steps: The first index 1 indicates always 
activation processes. If they are uncatalyzed, it is followed 
by u, if they are catalyzed by c; if they are effected by the 
solvent s, and if affected by a polymer p is added. Hence 
we get the classes (1u), (1c), (1s), and (1p). If there are 
several cases in one class, a third index 1, 2, --+ is added. 

10M. Polyani, Atomic Reactions (University of London 
Press, London, 1932); F O. Rice and K. K. Rice, The 
ror ceed Free Radicals (Johns Hopkins Press, Baltimore, 
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places as 

m,+c*—cm,* (1), 
or as mi+c —*m,*+c (2), (1c) 
and m,+cc—cm,*+c* (3). 


(1c, 1) or (1c, 3) seem to be characteristic for 
peroxide and radical type catalysts, while 
(1c, 2) appears to represent the action of metal- 
halide catalysts such as BF;,! AICI;, and 
SnCl,.t¢ 
- Collision with a solvent molecule has also 
been observed to produce nucleus formation, if 
the energy involved is sufficient. Again we have 
to consider two typical cases: 


The solvent mole- 
cule splits into 
radicals 


(1) 


(1s) 
The solvent polar- 


izes the monomer 
and facilitates 
opening of the 
double bond 


Finally we want to include activation by a 
vigorous collision with a terminated polymer m;, 
which might play a role in the later stages of the 
reaction, when much polymer is formed and 
which we can write as 


(1p) 


Altogether, we can summarize the different 
activation processes as follows: 


m,+s—*m,*+s. (2) 


(1u, 1)** 
*m*, (1u, 2) 
m+c*— cmy\*, (1c, 1) 
m+c (1c, 2) 
m,+cc >cm,*+c*, (1c, 3) 
m+s —*m,*+s, (1s, 1) 
(1s, 2) 
mi +m ;>*m,* (1p, 1) 


1 R. M. Thomas, W. J. Sparks, P. K. Frolich, M. Otto, 
and M. Muller-Cunradi, J. Am. Chem. Soc. 62, 276 (1940). 

122G. Williams, J. Chem. Soc. London 1938, 246, 1046 
(1938) ; 1940, 775 (1940). 

t It has to be mentioned that (1c, 1) and (1c, 3) are not 
true catalytical processes, because the catalyst takes part 
permanently in the reaction, and—as found experimentally 
—appears at one or both ends of the final chain-molecules. 
(1c, 2), however, might be regarded as catalysis according 
to the usual definition of the word. 

** According to this nomenclature, we shall use the 
following notation for the rate constants: Reaction (1, 1) 
_ the rate constant '&),, reaction (1s, 2) the constant 

1s, etc. 
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It can be seen that two altogether different 
types of nuclei can be formed. In case (1c, 1), 
(1c, 3), and (1s, 2), we get a particle having an 
unpaired or odd electron, which can only be 
disposed of by pairing up (saturating its spin) 
with another odd electron provided by another 
particle (atom or radical). In all other cases, the 
activated nucleus contains two unpaired elec- 
trons, which have been brought into this ex- 
cited (triplet) state by transfer of sufficient 
energy during the collision. Such activated 
nuclei (biradicals) are much more accessible to 
deactivation than the radicals (1c, 1), (1c, 3), 
and (1s, 2). Any collision which removes suffi- 
cient energy may bring the molecule back in its 
singlet ground state (closed double bond). It is 
true that the @ priori probability of such a 
triplet-singlet transition is small, but in the 
presence of perturbation which is always ma- 
terialized during a collision, the frequency factor 
may well be brought up from a value of around 
105 to its normal value of 10" or 10". 

It is therefore to be expected that nuclei of this 
type are much more sensitive to deactivation as 
compared with the radicals (1c, 1), (1c, 3), and 
(1s, 2), and it seems advisable to distinguish 
between deactivation by energy removal and 
termination by electron pairing. It will be shown 
‘later that experimental evidence favors this 
distinction and that accordingly we have to 
consider two types of inhibitors: those which 
catalyze the singlet-triplet transition, and those 
which take care of an unpaired electron. 

Finally it should be noticed that (1u, 2) will 
also be comparatively difficult to deactivate, 
because the closure of the two free valencies 
involves either the migration of a hydrogen atom 
or the formation of a four membered ring. 
Both processes require additional energy supply 
and are therefore presumably less frequent than 
the simple closure of the activated double bond 
in (1u,1; 1c, 1; 1s, 1; and 19, 1). 

We shall, therefore, introduce the following 
deactivation processes* which may destroy ac- 
tive centers before they even have a chance to 
grow into a short chain (such as a dimer or 


* As the first index 2 is already in use for propagation 
reactions, we suggest maintaining this use to characterize 
termination processes with 3 and deactivation processes 
with 4. The symbols u, c, s, and p are used as above. 
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trimer) : 
(u) (Rau) 
*mi*+c —m+¢, (c) (Rac) 
(4) 
*m,*+s —m+s (s) (Ras) 
+m (p) (Rap) 


In these processes pairs of asterisks attached to 
one molecule are destroyed, which corresponds 
to the closure of the double bond. 

Propagation in the form of chain growth, takes 
place by monomer addition to any of the active 
centers, as indicated by an asterisk in such a way 
that the number of asterisks on both sides of the 
equations are equal: activation is maintained 
during chain growth. 

Such processes are: 


* * (g, 1) (Ray) 
cm 541 (2) 
+m, 
2k 


(2g, 1) represents growth of the chain in one 
direction in such a way that a fragment of the 
catalyst appears at one end of the chain, while 
(2g, 2) describes chain growth in two directions. 
It seems that both types have been observed 
experimentally." 

There are, however, other, less simple processes 
which also leave the number of asterisks on the 
two sides of the equations unchanged and hence 
can be classified as propagation steps. They are: 

(a) Chain transfer. 


(t, 1) (*Ree) 
(2) 
(t, 2) 


It has been found experimentally that chain 
transfer can be catalyzed by certain substances 
(i.e., CCl,) in such a way that one part (¢) of 
the transfer catalyst (tr) terminates one chain, 
while the rest of it (r) originates another. This 


cm 


*m 


13]. Abere, G. Goldfinger, H. Naidus, and H. Mark, 
Paper presented at the American Chemical Society Meet- 
ing in Buffalo on September 9, 1942; H. Pfann, Thesis, 
Brooklyn Polytechnic Institute (May, 1942); C. C. Price 
and R. W. Kell, J. Am. Chem. Soc. 63, 2798 (1941); 
C. C. Price, R. W. Kell, and E. Krebs, J. Am. Chem. Soc. 
64, 1103 (1942); C. C. Price and D. A. Durham, J. Am. 
Chem. Soc. 64, 2508 (1942). 


KINETICS OF POLYMERIZATION 383 


can be represented by: 


cm *+trcm t+r*, (tr, 1) (Retr) 
(2) 
*m *+tr*m t+r*. (tr, 2) 


(b) Branching. As we do not have more than 
preliminary experimental indication for branch- 
ing of polyvinyl derivatives and no quantitative 
measurements as to its amount, we shall not 
include branching in this publication, but wait 
until further experimental development allows 
us to check formulas which can be derived in a 
similar way as shown for polymerization without 
branching in Part III. 

Termination or cessation destroys single active 
centers by providing for their odd electrons to 
pair up with other impaired electrons. The most 
frequent elementary processes of this type are: 


*m t+ (u) 
cm;*+c*—cm;c, (3) (c,1) 
cm +om,.* 0m 5440. (c,2) 


In these processes, pairs of asterisks disappear, 
but in such a way that each of them is con- 
tributed by another particle. This corresponds 
to the pairing of odd electrons belonging to 
different molecules. In case (3u), two chains 
saturate each other with one of their ends, hence 
two asterisks are still free on the right-hand side. 
If all four active centers would disappear, we 
would obtain mj;,., which represents a saturated 
ring having monomers. 

Termination can also be effected by certain 
molecules 6, which have the capability to dis- 
pose of an unpaired electron, without providing 
another unpaired electron (aromatic hydroxyl 
and nitro compounds, etc.). 


cm ;*+b—cm ;+, (b, 1) (*ks») 
(3) 
*m *+b—m +0. (6, 2) (hss) 


There seem to be still other cessation reactions 
possible, such as exchange of hydrogen atoms or 
reaction with the wall or with suspended par- 
ticles, but we shall not include them in the follow- 
ing calculations, as there is little or no quantita- 
tive experimental evidence for them at present. 


Ill. SOME KINETIC CONSIDERATIONS 


Activation in general will be given by 


dcm,* 
+k (Sa) 
dt 
d*m,* 
di =m, kiss +'Ripp), (5b) 
d*m,* 
(Sc) 
dt 


where, on the right-hand side the different 
contributing terms are enumerated. If one, or 
two, of them are preponderant, it will be justi- 
fied to neglect the others. As soon as an active 
nucleus is produced, two processes (in general) 
compete for it: the deactivation, which destroys 
it again, and the propagation, which builds it up 
into a chain. No deactivation,* however, can 
affect the nuclei produced by the terms ('h;.), 
(*k1-), and (*k;,) in (1c) and hence, we have to 
express the rate of deactivation by 


d*m,* 


dt 


= *my* +Rapp) (6) 


Propagation is represented by 


* 
dcm j+1 
dt 
d *m 


dt 


=k cm;*mi, (7a) 


(7b) 


and chain transfer by 


(8a) 


*m ;* (8b) 


It will be noticed that the rate constants 
(ke) and (?ke,) for normal propagation have 
been assumed to be equal, namely (k2,). This 
means that propagation is assumed to be inde- 
pendent of whether the chain grows in one or two 


* Note that deactivation involves the disappearance of a 
pair of asterisks from the same molecule. 
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directions and also independent of the chain 


length. There seems to be experimental support 
for this simplification. The same is true for 
chain transfer. 

Termination finally can be expressed by 


dcm;* 
(9a) 
t 
d*m;* 
t 


From these equations, we obtain for the 
steady state concentration of *m,* (active nuclei 
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being sensitive against deactivation) 
my c+ kiss +1 Ripp) 
1 => 
my (Rog+Rau) +Racc+Rass +Rapp 


and for cm,* (active nuclei being unsensitive 
against deactivation) 


(10a) 


Roy 

Furthermore, we get for the steady state 

concentration of active centers, incorporated in 


growing chains having already attained a certain 
length 


cm,* = (10b) 


(Rimi + Ries +R in) 


*m*=*m* = (11a) 
(Rau *m* (Romi +Rapp) 
lp ct +2k CC 
spb) 
Solving these two equations for *m* and cm*, respectively, one obtains 
1 22 2 t+ kict hist \3 
1 
cm* +4 |! — }. (12b) 


In the following we shall abbreviate the bracket 
term in (12a) with B(a) and the bracket in (120) 
with B(d). 

From these steady state concentrations of the 
growing centers we can now derive two experi- 
mentally measurable quantities, namely, the 
initial rate R with which the monomer is con- 
sumed (over-all rate of the reaction) and the 
initial number average degree of polymerization 
P., of the polymer formed. 

Considering that the rate of monomer con- 
sumption is proportional to the product of the 
concentrations of the monomer, and of the active 
uncatalyzed centers we obtain for the un- 
catalyzed reaction: 


R(initial, uncatalyzed) 


k 
(13a) 


3u 


and for the catalyzed reaction : 
R(initial, catalyzed) 


Rog 


=m (13b) 


3c 


The number average degrees of polymerization 
are 


P,,(uncatalyzed) 
k 


3u 


= , (14a) 


Rog 


+°k1.cc+——B(b) 


3e 


P,,(catalyzed) = . (14b) 


2c 


= 


Ja) 


ive 


a) 


b) 


a) 


b) 


D) 


2) 
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IV. DISCUSSION OF SOME 
EXPERIMENTAL RESULTS 


It is not the purpose of this paper to discuss 
fully the techniques and results obtained in our 
experimental investigations. We wish merely to 
indicate the general agreement between some of 
the more representative data and the theory as 
outlined above. 

Figure 1 shows the intrinsic viscosities of 
initially formed polymers for three initial mono- 
mer concentrations in toluene at 100°C plotted 
against the reciprocal square root of the catalyst 
concentration as added. It can be seen that the 
four measured points for the lowest (10 percent) 
monomer concentration are arranged with rea- 
sonable accuracy on a straight line, while for the 
two higher concentrations (20 percent and 40 
percent) the straight line bends toward the 
abscissae at lower catalyst concentrations. 

It seems that the equations derived above 
permit the following interpretation of this 
behavior : 

Assuming first that the intrinsic viscosity 
(over the limited range under consideration) is 
proportional to the weight (viscosity) average 
degree of polymerization, which in turn is 
proportional to the number average, we expand 
equation (146) and neglecting terms of higher 
order obtain 


kogm+'k3-c* 


P,,(catalyzed) = (15) 


Assuming further that the concentration of 
retarders or inhibitors of the type b is small we 
can further simplify (15) and get 


’ (16) 
+'k3-c* 


P,(catalyzed) = 


At very high catalyst concentration [very low 
values of 1/(c,)#] this reduces to unity, which 
corresponds to a very small value for the in- 
trinsic viscosity. As a matter of fact all three 
curves of Fig. 1 intersect the ordinate at values 
below »;=0.04. At lower concentrations of the 
catalyst the shape of (15) depends upon the 
magnitude of all three quantities on the right 
side. In general ke, will be large compared with 
‘kor, otherwise no high polymer would be formed 


at all. At small monomer concentrations we can 
therefore neglect 'k2,m; but not k2gm, at the side 
of 'k3,c*. Then the relation 


Rogm, 


P,(catalyzed) =1+ 


(17) 


is obtained, which calls for a straight line in the 
plot of Fig. 1, where c* is proportional to the 
square root of the catalyst as added. It can be 
seen that, in fact, the curve for the lowest mono- 
mer concentration employed is nearest to a 
straight line. At higher monomer concentrations 
'keym, cannot be neglected and chain transfer 
makes itself felt. It reduces the degree of poly- 
merization and with it the intrinsic viscosity of 
the polymer initially formed below the value as 
extrapolated from the initial straight lines. The 
curvature of the curve in Fig. 1 permits the 
determination of the rate constant for chain 
transfer, if the constants for propagation and 
termination are known. 

At very high catalyst concentrations the 
concentration of the monomer should have no 
influence on P, and all three extrapolated 
straight lines should intersect the ordinate at the 
same point. This is not exactly true in the graph, 
but is at least indicated with reasonable ap- 
proximation. 

Figure 2 contains the initial rates of polymer 
formation (which is equal to monomer consump- 
tion) as functions of monomer concentration at 
three initial catalyst concentrations in toluene 
at 100°C. Equation (135) shows that omitting 
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terms with 3, this initial rate should be 
R(initial, catalyzed) 


Rogm, 

2 

At low monomer concentrations the second 
term under the square root can be neglected and 
the over-all rate reduces to ki,mc*, while at 
higher values for m,, the first term can be neg- 
lected and one obtains 


RGinitial) myic**, (19) 
8c 


Plotting R/m, against m;} one is expected to 
have a straight line, the deviations of which 
are caused by the neglections in going from 
Eq. (136) to Eq. (19). Figure 2, being such a 
graph, shows that, in fact, at higher monomer 
concentrations (20, 40, and 50 mole percent 
monomer) a straight line is reasonably ap- 
proximated. At low values for mj, however, 
noticeable deviations take place. 
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Figure 3, finally, shows the influence of 
catalyst concentration. Plotting R divided by 
(c*)* against (c*)! one should again expect a 
straight line from Eq. (19). In Fig. 3, R 
divided by the fourth root of the catalyst as 
added is plotted on the ordinate versus the fourth 
root on the abscissae for two relatively high 
monomer concentrations (20 percent and 40 
percent). The points fall with fair approximation 
on straight lines. 

This approximate comparison of our experi- 
mental values with Eq. (14) and (19) together 
with the previous results of other authors as 
referred to in footnote 2 seem to indicate that 
the steady state approximation can be legiti- 
mately used to describe the chain mechanism of 
styrene polymerization if the various parameters 
(monomer concentration, catalyst concentration, 
temperature) are varied in limits as indicated 
above. 
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The (isothermal) Young’s modulus of ideal rubber is determined by statistical methods. 
The macro-state of a chain-molecule is here defined in terms of its length and not in terms of the 
distance between the ends of the chain. Section 1 deals with thermodynamical results holding 
for ideal rubber and in Section 2 the treatment based on ends-separation of the molecule-chain 
is described. Section 3 treats the problem on the basis of the length of the molecule-chain. 


1. 


HE marked difference between the elastic 
properties of rubber and a metal is well 
known, for not only are the elastic moduli in the 
case of rubber of a far smaller magnitude than 
those for a metal, but also, whereas in the case of 
a metal the tension of a stretched wire decreases 
with rise of temperature, in the case of rubber 
when heated in the expanded state the tension in 
general increases. The behavior of rubber is 
similar to that of a compressed gas for which the 
pressure (at constant volume) increases with 
temperature. This similarity between a gas and 
rubber is not merely superficial, but is the ex- 
pression of the underlying fact that in both cases 
the elastic behavior is primarily due to the 
thermal motions of the molecules. For the pur- 
poses of a statistical treatment a molecule of 
rubber is considered as a chain of a large number 
of links—the links being free to rotate about their 
junctions. Because of the thermal motions any 
molecular chain will be in constant sinuous 
motion; the shape of the chain will be changing 
from moment to moment just as the position of a 
molecule in a gas because of-the thermal motion 
changes from moment to moment. 

The basis of any statistical treatment is to 
enumerate the number of complexions corre- 
sponding to a given macroscopic state (i.e., the 
number of micro-states which all generate the 
same given macro-state), and then to derive 
entropy with the help of Boltzmann’s relation 
connecting entropy and probability. Several in- 
vestigators particularly Mark,’ Pelzer,? and 


1H. Mark, High Polymers, Interscience Publishers, New 
York, 1940), p. 69. 

2H. Pelzer, Reports on Progress in Physics (Physical 
Society, London, 1940) Vol. 6, p. 330. 


Wall* have made important contributions to this 
subject. In their treatment the macro-state of the 
molecule is defined by fixing the distance between 
the ends of the molecule-chain. In this paper we 
shall follow a different procedure in that the 
macro-state 1s defined not in terms of the ends- 
separation but, what appears more reasonable, in 
terms of the length (or extent) of the chain. In 
Section 2 the treatment (according to Wall) based 
on ends-separation is summarized, and Section 3 
deals with the treatment based on the length of 
the chain. 

In theoretical investigations on the elasticity of 
rubber it is usual to consider an ideal rubber 
which is defined by the property that 


(aU/al)r=0, (1) 


where U is the internal energy and / the length of 
the rubber piece. This relation corresponds to 
Joule’s law for the case of a perfect gas which 
states that the internal energy of an ideal gas is 
independent of the volume and depends only on 
the temperature. The approximate validity of 
Joule’s law for a rubber sample has been es- 
tablished by the investigations of Meyer and 
Ferri.* From the thermodynamic point of view it 
is significant that gases and rubber both approxi- 
mately obey Joule’s law. An immediate conse- 
quence of this is that if a rubber piece of 
unstretched length J» be stretched to a length / by 
subjecting it to a tensile stress I, then, [ and / 


3F. T. Wall, J. Chem. Phys. 10, 132 (1942); also E. 
Guth and H. M. James, Ind. Eng. Chem. 33, 624 (1941); 
34, 1365 (1942); Phys. Rev. 59, 111 (1941). 

4K. H. Meyer and C. Ferri, Helv. Chim. Acta 18, 570 
(1935); also R. L. Anthony, R. H. Caston, E. Guth, and 
L. E. Peterson, Ind. Eng. Chem. 34, 1349 (1942); J. Phys. 
Chem. 46, 826 (1942). 
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must be connected by a relation of the form 
r=Tf()), (2) 


where T is the absolute temperature and f(/) is 
some function of length / only. For, if dS denotes 
the increase in entropy, we have 


dU=TdS+Tdl, 


as 


and using (1) 


or 


If we use the thermodynamic relation 
as or 
OT / ; 


r=Tf(l). 


As T depends only on the change in length, we 
may write 


we have* 


P=Tf) =TF(l/lo), 


where J) is the unstretched length. When the 
stress is small T will be proportional to (1—1,)/lo 
and, thus in the region where Hooke’s law is 
valid, (2) reduces to 
T=T7R’ 


where R’ is a constant. The object of the sta- 
tistical theory is to determine: 

(i) the value of the constant R’ in terms of 
chain-length and the molecular weight of rubber, 
and (ii) the form of F(J/lo), particularly when 
1/lp is large. 


The starting point of the usual statistical 


theory of rubber is the probability relation (well 


* Thus, for ideal rubber the temperature-coefficient of 
increase of tension for constant length is 

** The length J) may be supposed to correspond to the 
term ‘‘d” in van der Waals’ equation for a gas. 


known in the discussion of Brownian motion) 
(x)dx=A ( (6) 
x)dx=A exp | ——— }dx, 
2rd 


where p(x)dx denotes the probability that the 
distance between the two ends of a chain of v 
links, each of length b=X/v, lies between the 
limits x and x+dx. The normalization constant 
A is given by 


2 


We shall consider throughout in this paper only 
the one-dimensional case, i.e., we shall assume 
the chain to lie entirely along the x direction. 


The generalization to three dimensions is easily — 


made, but that will not be discussed here. The 
object of the present paper is to describe the 
statistical theory of rubber in terms of the length 
of the molecule-chain, and to compare this ap- 
proach with the usual one which is based not on 
the chain-length but on the ends-separation. 
(The chain-length p should not be confused with 
\=vb which denotes the length of a fully ex- 
tended chain. Because of the folding of the chain, 
p in general will be extremely small compared 
to X.) 

Let us now consider an assembly of a large 
number N of similar molecule-chains. Let N(x)dx 
denote those chains of the assembly which have 
their ends-separation lying in the interval 
(x, x-+dx). To be precise, we cannot take x to be 
a continuous variable, for it will vary discon- 
tinuously in steps of b, and, therefore, to discuss 
the statistical distribution we suppose the inter- 
val of length \ divided into a large number of 
zones or cells such that the size of any cell is 
large compared to b and yet very small compared 


to X. Let the rth cell correspond to the interval x, 


to x,+dx,, and let N, denote the number of 
molecule-chains, having their ends-separation 
lying within the above interval, then, we have 


N,=N, 


and 
«,N,=EN, (8) 


where Z is the average ends-separation. 
Consider a distribution defined by the numbers 
(Ni, Ne, ---, Nr, -++). The number of ways of 
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realizing this distribution will be given by 


W=N! TI (9) 


r N,! 


when ~, is the probability of a chain being con- 


tained in the rth cell; i.e., p, is the probability 
that the ends-separation lies within the interval 
(x,, x-+dx,). We have from (6) 


pr=A exp Jas p-=1. (10) 
2nd 


We shall first determine the value of N, for the 
most probable distribution, i.e., the distribution 
for which W is a maximum. By differentiating (9), 
or even otherwise, it follows that for the most 
probable distribution 


N,=Np,, (11) 


where a bar over a symbol denotes its value for 
the most probable distribution. The maximum 
value of W is given by, 


(pr)¥?r 
W=N! 


From (9) and (12) we obtain, after using Sterling’s 
theorem, 


(12) 


Dy Nr 1 
W=N*(2xrN)} 
r \N, (2xN,)} 
and 
(2xN)} 


and, therefore, introducing Boltzmann's theorem, 
we have for the decrease of entropy 5S, 


W 


i". 


Np, 
=), (N,+2) log 


(M44) log (13) 

N, N, 
where k is the Boltzmann constant. As W is the 
maximum value of W, log (W/W) <0. Equation 
(13) is an important one as it gives the decrease 
In entropy in passing from the most probable 


distribution (sometimes referred to as W-distri- 
bution) to any other distribution. As an illustra- 
tion consider a distribution which differs from the 
most probable distribution only as regards the 
population in the rth and sth cells: N,=N,+1, 
N.,=N,.-—1, then, we have 


— (N,+1) 1 


N,+1 
N, 1 


The application of Eq. (13) to the case of 
rubber is made in the following manner. The 
unstretched state of rubber is identified with the 
state of maximum probability. The entropy of 
ideal rubber in the region of Hooke’s law is given 
by, using (4) and (5), 


( (14) 


as 
(—) dl=y(T)—R’ 
al/ 21 


S=(T)+ J 


0 


where ¥(7) is a function of temperature only. 
The entropy is thus, for a given 7, maximum for 
l : 

The probability W for the stretched state is 
determined by assuming that when rubber is 
stretched from its original length J) to l=alo, 
then, the ends-separation of every molecule-chain 
is also increased in the same ratio.* For the 
stretched state the rth cell will contain chains 
that originally in the unstretched state had their 
ends-separation lying in the range x,/a to 
(x,+dx,)/a, so that after stretching from 1) to 
l=alo, the ends-separation lies in the interval x, 
to x, +dx,. The value N,, therefore, is equal to the 
number of chains which in the most probable 
distribution (i.e., unstretched state) have their 
ends-separation in the interval x,/a to (x,+dx,)/a, 
and hence from (10) and (11) we have, 


2\3 \ dx, 
v.=(—) exp (— ) 
2rba?/ 


*We of course ignore the extremely small number of 
molecules whose ends-separation already in the unstretched 
state is comparable to i. 
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Thus we have 


= r a ’ 

which using (15), and replacing the sum by an 
integral reduces to 


N(o®—1) 2 
(2d)! 


x? \x? dx 
x f exp (- —, 
0 2r\ba?/ a? a 


2N(a?—1) 
= N log f 2 exp (—2*)dz, 
0 


or we have 


0 


Therefore, we obtain 


as 1 /] 
r--7(—) --ver|--—| 
lL Ie? 


NkT 


(17) 


The Young’s modulus (y) for the isothermal case 
is given by 


oT 1 / 
lL I? 


and if do denotes the density of rubber in the 
unstretched state and M the molecular weight of 
a molecule-chain, we have 


2RTdo 
Io M 


’ (18)* 


where R is the gas constant. 


3. 


We shall now develop, following Section 2 with 
necessary alterations, the statistical treatment by 


* Equation (18) is identical with (19) of Wall (refer- 
ence 3) 


defining the macro-state of a molecule-chain in 
terms of its length. It is improbable that a 
molecule-chain will be found in a completely 
unfolded condition, but, on the other hand, it 
will be folded at many of its junctions, and, 
therefore, in general, the length of the chain will 
be much smaller than its uncoiled length A. The 
probability p(p)dp that the length of the chain 
lies between p and p+dp is given by the following 
expression (when the number of links v is large) 
which has been derived by Daniels® in a recent 
paper: 


( \d 8 ( ( 


(19) 
f 


b being the length of each link, and \=vb. Using 
the properties of #-functions, Eq. (19) may be 
transformed into the following form which is 
more convenient for small value of p: 


b(p)dp =— { — —1)?} 


p® 
(2s 
Xexp ( (20) 


2p" 


The asymptotic expressions for the two limiting 
cases are, for p—0 


P(p)dp= exp (21) 
2p? 


and for p> 


8 
dp= ——- }dp. 22 
p(p)dp ( ~~) p. (22) 


It may be remarked here that, whereas the most 
probable ends-separation of a molecule-chain is 
zero, the most probable value for the length of the 
chain is 1.346(\d). Moreover, almost all the 
chains have their lengths lying in the interval 
(0.5(Ab)#, 3.5(Ab)#). The average value of p as 
obtained from (19) is p=2(2Ab/m)!=1.596(bA)?, 
whereas the average value of the ends-separation 


5H. E. Daniels, Proc. Camb. Phil. Soc. 37, 244 (1941). 
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is 
|x| = p/2. 

Let us consider a statistical assembly of a large 

number of chains at temperature 7, and let N, 

denote the number of chains having their lengths 

lying within the interval p, to p,+dp,. The 

number (W) of ways in which a distribution 


(Ni, No, ---, of the assembly can be 
realized is given by . 
=N! 
W NI (23) 
N-=N. (24) 


In (23), p, denotes the probability that the 
length of a molecule-chain lies in the interval 
(pr, pr+dp,), i.e., we write p, for p(p,)dp,. From 
(23) and (24) it follows, as in Section 2, that the 
most probable distribution is given by 


N,=N,=Np;; (25) 


and if W denotes the corresponding value of W, 
we obtain for the decrease of entropy S, 
S-8 N, (26) 
k r N,; 

When rubber is stretched from its original 
length 1) to l= alo, let us assume that the length of 
each molecule-chain is also increased a times.* In 
the stretched state the rth cell will contain chains 
which in the unstretched state, had their lengths 
between p,/a and (p,+dp,)/a. 

Therefore, we have 


r\@pr 
<4 (27) 
a 


a 


8a’ Ta 6 


a‘ 


and using (25), and (26) we obtain 


S-S ap(p) 
Je (‘) log 


ap(ap) 
=N (p) log ———-d 


log asf P(p)dp 


p(ap) 
+ ie 
J log 


= N log atN f b(p) log ep (28) 


Thus we have 


0 
as 
al/ 


NkT {1 ‘(ap) 
=— {-+ f pb(o)— do}, 


p(ap) 
or 
RTdoj1i 1 
(29) 
M \a a 
where J denotes the integral. 
b'(p) 
pp(p/a) 
0 P(p) 

p\ b'(p) p\?'(p) 
0 al p(p) p(p) 
(30) 


For evaluating the first integral we use the ex- 
pression for p(p) given by (20), and for the second 
integral we use that given by (19), and we obtain 


8a*(a®?—1) 3203 
2 2 


| (176 — 367+ 182")a?— 182 — 16 (724+32)a*— 36 144a? 
x 
(8+)? 


32a3(at—5a?+2) 8a? 


re! 


3r—-5 
(8+a?) +~exp (- + 


_2(a?— 1) | 
(8+a*)® 


9a? 27 
8(9a?— 1) a3 1) 
exp (402) 2 


2 


(31) 


*In making the above assumption, we ignore those molecule-chains which in the unstretched state of rubber 
already possess a length comparable to \. As almost all the chains have their length lying between 0.5 (Ab)+ and 
3.5 (Ab)!, the maximum permissible value of @ is of order (A/b)}. 
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where 
¢(x) f —dw, 
w 
and 
Erfe (x) = f exp (—w?)dw. 


In the case of a~1 we have directly from (29) 


l—ly 
lr =12.6-——_ ——, (33) 
M 
and for a>1 
RTd, | 
['~4 log 2 (34) 
d 0 
‘A 
12 
10+ 
RTdg 
B 
4 be 
2b 
t/t,» 
( 2 4 


Fic. 1. The curves in the figure show the relation between 
Young’s modulus y and the stretched length /. The curve 


A is based on chain-length, and the curve B on oni 


separation of the chain. 


Thus we have for the isothermal Young’s 
modulus 


or or 


RTdy 
12.6—— fora~1 
=) (35) 
RTdy 1 
4 log for a>. 


0 


When a<<1, i.e., for the case of high compression, 
Eq. (29) gives © 


—, (36) 
‘2M F 


y= 


The curve A in Fig. 1 exhibits the relation 
between y and as determined from Eggs. 
(29), (31), and (32). The curve B represents 
the same relation as given by Eq. (17) which 
is derived on the basis of ends-separation 
of a molecule-chain. The extension to three 
dimensions and the comparison of the theoretical 
results with experiment will be taken elsewhere. 
The purpose of the present paper is to show how 
the results of the usual theory are modified when 
the chain-length and not the distance between the 
ends of the chain is taken as the parameter to 
define a macro-state of the molecule-chain. 
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Gel Formation in Vinyl-Divinyl Copolymers* 


W. H. STOCKMAYER AND HOMER JACOBSON 
Department of Chemistry, Columbia University, New York, New York 
July 12, 1943 


HE gelation of polymerizing systems containing 
multifunctional reactants has been shown by Flory! 
to be due to the sudden formation of macroscopic network 
molecules. For polycondensations, in which the reaction 
mechanism is very simple, a statistical theory can be 
formulated which permits a rather accurate prediction! of 
the extent of reaction at which gelation occurs. In the case 
of vinyl polymerizations, on the other hand, a completely 
a priori prediction of the gel point is not feasible because 
the mechanism is complex and but imperfectly understood. 
Conversely, however, experimental observations of the gel 
point may be used to study certain features of the reaction 
mechanism, as we shall indicate below. 

In a polymerizing mixture of monovinyl and divinyl 
reactants (subject to the experimentally realizable condi- 
tion that all double bonds are equally reactive and hence 
enter the chains at random), we find that the expected 
number of cross-links to which a given cross-link is con- 
nected by a single chain is pp(Aw— 1), where p is the fraction 
of reacted double bonds, p is the over-all fraction of double 
bonds residing on divinyl units, and A. is the weight average 
number of reacted double bonds in a stable polymer chain. 
As Flory! has shown, as soon as this expectance exceeds 
unity there is a finite probability that very large cross- 
linked molecules shall form, so that the extent of reaction at 
the gel point is 

pe=1/p(rv—1). (1) 


Actually Eq. (1) gives a lower limit for p., since intramolec- 
ular reactions leading to cyclic structures are not allowed in 
the derivation, but the error is quite small. 

It is to be emphasized that in deriving Eq. (1) no as- 
sumptions are made concerning the reaction mechanism 
other than the aforementioned one that all double bonds 
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are equally reactive. This condition can be realized by a 
proper choice of reactants if, as has been demonstrated for 
polyesterifications? and is very likely in vinyl polymeriza- 
tions,* the rates of all steps in the polymerization process 
are independent of the size of the molecules on which the 
reacting functional groups are located. This state of affairs 
does not, as was erroneously implied in another paper,‘ 
necessarily lead to the same molecular size distribution 
found in polycondensation reactions, which involve only a 
single reiterated step. 

We note further that ‘branching’ reactions (e.g., 
transfer of a hydrogen atom from one free radical to 
another, resulting in a stable molecule and a biradical) are 
not excluded from the mechanism, so that the “chains’’ 
whose weight average length \,, enters Eq. (1) are not 
necessarily linear. They may be pictured as the aggregates 
that would remain if all the divinyl units were suddenly to 
be chopped in half. 

Since p, and p are easily measured, Eq. (1) appears to 
offer a new and very simple method of studying the varia- 
tion of \~ with temperature, concentrations, and extent 
of reaction. It has special interest for systems in which the 
rates are independent of p, for in such cases Aw is identical 
with the weight average degree of polymerization of the 
pure monovinyl system. One might presume such a condi- 
tion to exist if the divinyl units are essentially dimers of the 
monovinyl units (for example, vinyl acetate and divinyl 
succinate, methyl methacrylate and glycol dimethacrylate, 
etc.), since the effects of changing medium would then be at 
a minimum. Alternative ways of measuring \~ are either 
tedious (fractionation, sedimentation) or of unknown 
validity (viscosity). 

Equation (1) also applies, with a redefinition of symbols, 
to experiments in which previously prepared linear 
polymers are cross-linked by an added reagent (e.g., the 
vulcanization of rubber). In this case p signifies the fraction 
of units (in the polymer characterized by \) susceptible to 
the cross-linking reaction, and » is the fraction (at the gel 
point) of such units involved in cross-links. If the chains are 
long and of uniform length, Eq. (1) reduces to a result 
previously given by Flory.® 

We intend at a later time to present a more complete 
treatment of the theory, and to report on experiments now 
in progress. 

* Publication assisted by the Ernest Kempton Adams Fund for 
Physical Research of Columbia University. 

1P, J. Flory, J. Am. Chem. Soc. 63, 3083 (1941). 

2P. J. Flory, J. Am. Chem. Soc. 61, 3334 (1939). 
as 5 “4 . G. Herington and A. Robertson, Trans. Faraday Soc. 38, 490 


4W. H. Stockmayer, J. Chem. Phys. 11, 45 (1943). 
5P,. J. Flory, J. Am. Chem. Soc. 63, 3096 (1941). 
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